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SYNOPSIS 

SUBHASH CHANDRA RASTOGI 
M.Tech,(Mech.) 

Indian Institute of Technology Kanpur 
October, 1S73 

A 1HTH0D0L0GY FOR THE DESIGN OF AN ASSEMBLY LINE 
WITH VARIABLE TASK TIMES 

In the present study a design methodology which 
accounts for the following important features of 
assembly lino production has been developed : 

1, A case where the time requirement fpr at 
least one task is more than the cycle time, 

2, Variability in the task times, 

3, Failure characteristics of the work stations. 

The design methodology has been divided into 

two phases. The first phase considers only the first 
feature whereas the other two features are taken care 
of in the second phase. 

The first phase consists of allocation of the 
total number of the tasks of the line among the various 
work stations, in order to achieve maximum compactness 
and uniformity in the distribution of total work load. 
The task times have been considered deterministic in 
this phase. The case, where the cycle time is less 
than the task time for at least one task of the assembly 
line, has been incorporated in this phase by intro- 
ducing parallel stations. COMSOAL technique has be'en 
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used to achieve the maximum compactness and the 
concept of smoothness Index is used to achieve maximum 
uniformity* 

Tho second phase accounts for the variability 
in the task times as well as the failure of the work 
stations. Buffer inventories are introduced to account 
for both of those features. In this phase the optimum 
levels of the buffer capacities are determined which 
minimize the total cost of production per unit produced 
from the line . Digital Simulation and Univariate 
search teclmiqu.es have been used to achieve the optimum 
size of tho buffer capacities. 

The proposed methodology has been computerized 
and validated with the help of the data collected from 
an ei.ectronics industry manufacturing television sets. 



C HAPTER I 


INTRODUCTION 

■Jith the advancement of civilization and consequent 
enhancement in the standard of living of the people in 
general, there have been great changes in the demand patterns 
and functional requirements of all sorts of commodities* 

The demand curves for both civilian goods and military 
equipments have been rising in all its multi cudes* This 
posed a big challenge to the managers and technicians all 
over the world to revise the production system radically as 
to make them capable of meeting the increased demands of 
varied nature of commodities. Research scholars put tncxr 
heads together to devise ways and moans for achieving mass 
scale production. These efforts resulted in technologically 
Improved system suited for mass production* Line system of 
production is one such advancement whose value was recognized 
during the period of transition from small scale to large 
scale production. 

Henry Ford is credited for being the first person 
to implement the line system of production in its highly 
developed form for Ford Highland Park Plant. This was in 
ISIS • It was observed that this system brings out a 
sharp increase in productivity and also improves the 
quality of the product at reduced cost of production. 
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The method was found so attractive and useful that it has 
been adopted by all kinds of industries ranging from toys 
to automob ile s . 

Tie line system of production, as Demyanyuk £s]put s 
it, is a production system where the machining or assembling 
operations are allotted to definite pieces of equipments and 
operators. The operations required by the product are 
arranged according to their precedence only. The semi- 
finished products are moved from one work station to the 
next after completion of the preceding operation. 

Among the problems associated with line system of 
production, the line balancing is the most important one. 

The essence of line balancing problem lies in apportioning 
the total work load required for assembly work among the 
work stations uniformly and compactly following the 
precedence restrictions. 

This problem generally involves either, 

(1) minimization of number of stations for a 
specified output rate, 

or (2) maximization of the output rate for a fixed 
number of work stations* 

A number of analytical and heuristic models have 
been suggested for the purpose of balancing a line for a 
required output rate. A few research workers have concerned 
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theinselve s to the second criterion* 'the models available 
in literature are signif.ica.ntly different from the real life 
situation of an assembly line, l-'oreovcr, the objective 
function in most of the ava.il able moo.cls is the minimiza - 
tion of total idle time of the assembly line. This does 
not guarantee the minimization of total cost of production, 
which is of greater interest to the industries* The factors 
which have been either overlooked or have not been paid due 
attention are listed as xollows t 

1. The random variability of the station times ■ 

2. Uncertainties related with the failure of the 
stations 

3. A situation where at least one task time is so 
large in its time content that it cannot be 
performed by a single operator to fulfil the 
required output rate. 

These factors very much govern the actual throughput 
and cost of operation of the line, and should be taken into 
account to represent an actual line realistically. The 
present work has been directed to develop a model which 
considers each of the above listed factors. Interstage 
inventories have been introduced in between ea.ch pair of 
successive stations to absorb the effect of variability of 
the task times. The problem has been divided into two 
phases : 
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Phase I : Maximization of output rate by balancing the 
line, ignoring the effect of variability of 
task tine. 

P hase II i Do termination of optima], levels of interstage 

inventories in between every pair of successive 
stations, taking into account the variabilit}? - of 
task tine and failure of stations* 

The model, incorporating above listed two phases, 
has been programed in FORTRAN IV language for use on IBM 7044 
digital computers. The model has been applied to an assembly 
line of J.IC* Electronics, Kanpur. Although the model has 
been primarily tested for the said television ind.ustry, 
every care has been taken to make the model general in 
nature and, therefore, may be adapted for other assembly 
lines without much difficulty. 



CHAPTER JJ. 


LITERATURE REVIEW 

The need for balancing the production line was 
felt as early as in 1900 but a successful technique to 
achieve this was not found until 1913, Since then much 
efforts have been put in this direction and as many as 
twenty models have come up during the last couple of 
decades. Some models tal:e a realistic approach to solve 
practical problems, while others deal with its theoretical 
structure with a view to gain deeper understanding of the 
problem without making any attempt towards their applica- 
tion to i arge scalo balancing problems. 

The first paper published in this field was by 
Salve son (s3j in 1955, though credit should have g-one to 
Bryton |4] , who was first to analyse this problem 
analytically in 1954. Bryton 1 s research work was not publi- 
shed however. Sa-lveson’s paper suggests that one should 
enumerate all possible work stations and then use Linear 
Programming to choose the combination of stations that 
best meets the requirement placed on the solution. The 
problem was again formulated as a Linear Programming 
Problem (this time as an Integer Programming Problem) 
by Bowman j^s] , in I960. Salve son* s method is unacceptable 
for it results in fractional value of work stations which 
cannot be achieved in reality. This difficulty is overcome 
in Bowman’ s model . 
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Bryton’ s methodology suggests interchanging the 
pair of tasks between the largest and the smallest stations 
(on time scale) according to a specified rule, to attain 
a better balance. The simple concept of each transfer 
causing improvement seems rather unconvincing, and he 
himself admits it. Jackson’s enumeratlve procedure [is] , 
published in 1G56, promises an optimal solution, but it 
requires enormous computational efforts while applying to 
real life problems. In fact, all the models mentioned 
above need too much of computational efforts to be practical 
and therefore, are of little value to industries. 

In 1961, Kilbridge and Wester [l9j proposed a 
technique which is found useful for a variety of line 
balancing problems. It is a heuristic technique and requires 
■little ingenuity and little knowledge of arithmatic# The 
method makes a systematic use of the precedence graph and 
appears to be a useful way for arriving at a balanced line 
by hand computations. Comparison of its ability with other 
techniques has been done later in this chapter. 

In the same year, another procedure, popularly 
known as RPW (Ranked Positional Weight Technique), was 
proposed by Hclgeson and Birnie. jli] . The method is 
logically very simple and needs only one iteration. It 
is workable for small lines but significantly deviates 
from the optimal solution for large size problems. An 
improvement to this method was suggested by Mansoor(24j in 
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1964, He claims that an optimal solution is guaranteed 
if proposed modifications are followed, but the amount of 
work required may he quite large. 

In 1963, Held, Karp and Sherishian jisj developed 
another ana.lyt.ical procedure using Dynamic Programming 
approach. This has a limitation on the size of the problem 
it can handle, o then-rise it is the best technique available 
so far, for line balancing problems. To cope up with the 
large size problems, the authors have suggested a method to 
group some tasks and consider them as a single task. This 
modification reduces the size of the problem, however, it 
does not guarantee an optimal solution. The capabilities 
of this technique are discussed in details later in this 
chapter. 

In 1963, another model, called Successive Maximum 
Elemental Time Technique, was introduced by Hoffman 
This model develops a somewhat modified form of precedence 
matrix which is then used to generate all feasible permuta- 
tions of tasks and assigns them to the various stations. 

It gives only an approximate solution, as was evaluated by 
Ignall [l7] and others [si, 27] . 

In 1965, Tonge [35] came out with some probabilistic 
combinations of heuristics for choosing the next task to be 
assigned to the station. Two other Important methodologies, 
for balancing the line, were developed by Moodie and 
Young [2S] , in 1965, and by Arcus[l] , in 1966. It is 
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hope 1 ? that the former technique provides good assembly 
line balances for both manual and computer mode of solu- 
tion, but not much pain has been taken by the research 
workers to verify this, while other methods have been 
tested for their effectiveness by more than one author* 

.arcus r s method, popularly known as COMSOAL (Computer Method 
of Sequencing Operations for Assembly Line) , is a heuristic 
method and is claimed to be best for large size line 
balancing problems {26| . 

Some other methods for line balancing were also 
suggested by Hu [l6] , Klein [22] and Gufahr and Nemhauser [12J , 
but were of little practical value because of innumerable 
computations required, 

Ignall [l7] , in 1965, and Master (27) , in 1970, 

conducted a comparison study of some of the prominent line 
balancing techniques. Their work is one of the most valuable 
contributions to the area of line balancing because no 
attempt had been made in the past to compare the various 
line balancing techniques. Master [27] considered the 
models proposed by : 

1. Helgeson and Birnie 

2. Hoffman 

3 . Arcus 

4. Held, Karp and Shershian 

5. Kilb ridge and Wester 
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Besides these, he also considered some of the 
commonly used heuristics, for generation of sequence of 
tasks* The heuristics, he considered for the purpose of 
line balancing are : 

1. Lexicographic Order 

2. Humber of immediate follower task rule 

3. Task time ordered rule 

4. Random sampling rule 

Raster’s investigations revealed that, the models 
of Held and Arcus consistently stand out as .the best, 
requiring the least Idle time for any number of stations. 
But Held 1 s technique seems to have limitation on to the 
size of the problem to which it is applicable , otherwise , 
as Raster claims, it is the best performer. Lexicographic 
Order Rule was found to be worst throughout. Master 
concluded that only COMSOAL by Arcus, and, heuristic line 
balancing technique by ICilb ridge and. Hester are the most 
effective techniques for large size practical problems 
(having more than 40 tasks). For moderate and small size 
problems, Held’s technique produces the best results. 

All the methodologies mentioned above do not 
consider the feature of variability • of task time, which 
is a very important factor in deciding the actual output 
from the line. The time requirement for each task Is 
generally not constant and may be represented by some 
amperically determined probability distribution. This 
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problem was taken up by some authors, namely, Hoodie and 
Young j2s] } Mansoor [ssj and Ians in g and Downing [32j . 

Hoodie and Young tackled the variability feature 
by alloting some allowances to the operators, so that 
actual station time does not exceed the designed cycle time, 
for a given confidence level. Provision of allowance to 
the operator reduces the output rate from the line and 
therefore, does not assure to obtain an optimal solution, 
so long as the objective function concerns with the cost 
aspect of the line. 

Ramsing and Downing [32j also tried to incorporate 
the variability feature of task times in the manner which 
was suggested by hoodie and Young. But their treatment is 
different from that of the latter in the following ways. 
Ramsing and Downing used sampling technique to determine 
the task times, for a given confidence level, and then, 
with the sampled values of the task times, they balanced the 
line using RFW technique. These authors also, like Hoodie 
and Young, did not consider the possibility of introduction 
of buffer inventories, as was later suggested by Freeman [9]. 
As per Freeman and others [ 2 , 3 ] introduction of buffer 
inventory is more economical way of achieving improved 
productivity. 

Mansoor [ssj , also introduced another model, 
which allots allowances to the operators. Along with this. 
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lie has also suggested, to introduce an incentive plan for 
the operators. It might he said that the incentive plans 
suggested by them, while convenient for their model, seems 
to he totally impracticable in real life setting. 

None of these authors [26, 28, 32j consider 
interstage inventories, paralleling of stations, failure 
of stations, etc. The fact is that all the models suggested 
so far, whether incorporating variability of task times or 
not, stress more upon reduction of total idle time instead 
of minimizing total operating costs. The present work is 
an attempt to develop a cost model of the assembly line, 
which also considers the following features of the line, 
to represent it realistically : 

1. Variability In the task times, 

2. Failure characteristics of the work stations and, 

3. A situation where, at least one task Is so 
large in Its time contents that it cannot be 
performed by a single operator to fulfil 

the requirement of minimum desired output rate. 



CriAPTER III 


PRC ELS II -DESCRIPTION AND FORMULATION 
3 .1 TERMINOLOGY 

Tiie following terminology is generally used in 
literature, in the context of line balancing problems [loj : 

1. Task : A feasible subdivision of work content of the 

total assembly. 

2. Task Time : Time content required to perform a task* 

3. Work Station : A combination of .tasks that must be 

performed serially forms a work station, 

4. Work Station Time : The time requirement for completion 

of total work load assigned to a work station. 

5. Cycle Time : The maximum amount of time available to an 

operator to produce one unit. 

6. Precedence Relationship : The description of the tech- 

nological order In which the various tasks must 
be performed in the assembly operation. 

7. Line Length : Total number of work stations on an 

assembly line. 

8. Blocking : A station is said to be blocked, 

a. if work on the item is completed but It 
cannot be passed to the next station 
either because the next operator is not 
free, or there is no available space for 
the inprocess inventory, and 
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b. If an operator, though free, Is unable 
to serve an Item either because of the 
exhausted inprocess inventory or the 
preceding operator has not finished the 
work assigned to him. 

9. Buffer Inventory : Semi-finished items stored between 

the two work stations to lessen blocking. It 
is also termed as inprocess inventory. 

10. Slack Time : This is the idle time at a particular 

station. It is defined as the difference of 
the cycle time and the work station time of a 
station, 

11. Balance Delay : It is a measure of the total Idle time 

on the line due to the imperfect division of 
work among stations • 

3.2 GENERAL DESCRIPTION 

The problem of line balancing is essentially to 
distribute the total work load of the assembly line among 
the operators as uniformly and compactly as possible, 
keeping in view the desired output rate and the precedence 
restrictions imposed by the various assembly operations to 
be performed on the product. 

The existing models till to-date, fail to account 
for the following aspects of an assembly line problem*. 
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Thev are : 

(1) An operator cannot follow, in general, the times 
specified by the design* The involved human element 
brings about an uncertainty in the work station times. 
Hence, the consideration of the work station time as 
deterministic value is inadequate to represent the 
system realistically. The work station time should 
be considered as a stochastic variable following a 
probability distribution function* 

(2) The variability in the time taken by an operator on 
a station results either in idle time of an operator 
or a high level of ihprocess inventory. Both of 
these result in a reduced output compared to the 
value for which the line has been designed. The 
problem can be solved by introduction of optimum 
buffer inventories to account for the variabilities 
in the task times. An infinite level of buffer 
inventory between each pair of stations will remove 
the blocking altogether hence increases the output 
rate to the maximum. But, it is neither economical 
nor feasible to have an infinite lovel of buffer 
inventory. On the other hand, a zero buffer inventory 
will cause maximum blocking and will reduce the 
output rate to its minimum. Therefore, the need for 
obtaining a compromise between the two extremes is 
felt so as to optimize the given objective function 
for the assembly line. 
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There are two factors, which govern the buffer 
capacity sizes of the line. They are ; 

1. Variance of work. station time, and 

2. Buffer holding cost. 

The optimization models, reported in the literature 
for the assembly line, have assumed equal variances of the 
station times and equal buffer capacities in between each 
adjacent pair of stations* In reality, the variances of 
the work station times need not be equal. The amount of 
variance at a station depends upon the nature of the work 
assigned to the work station. Therefore, it should not be 
taken as equal for all the stations. 

Most of the models developed so far, consider 
cost due to buffer inventory based on a buffer cost value 
averaged out over the whole line. In a real life situation, 
the value of the semifinished product increases as it 
passes through the line and therefore, the holding cost 
associated, should also be considered in accordance with 
the change in values of the product. This variation in 
values, and hence buffer holding costs, becomes more signi- 
ficant for long assembly lines, having sometimes the 
difference of values, between initial and 'final stage as 
high as ICO to 200 times* 

(3) The deterministic models of the line balancing also. fail 
to account for the failures of the work stations. The 
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failure of a work station is also a stochastic variable 
following a particular probability distribution. The 
failure of a work station very much affects the output 
rate of a system. The effect of a work station failure 
can be nullified by choosing an optimum level of 
buffer inventory. 

(4) Most of the models reported in the literature assume 
that the task times for all the tasks of the assembly 
are less than or equal to the cycle time. But, there 
are cases where there exists at least one task in the 
assembly, whose task time requirement cannot be met by 
a single work station within the specified cycle. A 
general model of line balancing should incorporate this 
feature. 

3 .3 THE PROBLEM 

Before delving into the problem formulation proper, 
it is essential to have a better understanding of the system 
under study, so as to help modelling. Therefore, a descri- 
ption of some of the important features of the system is 
presented here. 

There are m tasks in the whole assembly line. 

These tasks are to be assigned among the n work stations 
uniformly and compactly, following the precedence restric- 
tions. The assignment of tasks to the various work stations 
is done by the technique explained in Phase I of Chapter TV. 
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The task times, and hence the work station times are 
random variables, therefore, buffer inventories are intro- 
duced to absorb this variability. • 

A general assembly line can be represented schema- 
tically as shown in figure 1. 


Unlimited 
supply of 
items 





Pifl.-l 

In figure 1, : indicates the work station i and, 

: shows the buffer space between 

stations Sj_ and s i + i 

S n - in the figure 1 represents the work station i, loaded 
with the various tasks. 

The following assumptions have been made here : 
1. There . is an unlimited supply of items for the first 


work station. 

2, There is an unlimited storage space following the last 


work ' station. 

The first work station, S-j_, picks up an item 
from the supply point, performs the operations required on 
it and then passes it on to the next work station Sg« 
Similarly, the Item passes through every work station, in 
the prescribed sequence. Every work station adds value 
to the Item and then passes it to the next work station* 
The process is continued till the finished item comes out 
of the line. 
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Tho assembly lino, shoun in figure 1, can be 
considered as a single channel queuing system with n 
subsystems, each consisting of a service facility and the 
.items in the queue, waiting for the service. The work 
stations a. re service facilities subject to the random failure 
and the buffer inventories are the queues. There are two 
possible events which change tho state of tho subsystem : 

1. Arrival of an item in the subsystem, 

2. Departure of an. item from tho subsystem a f ter 
being serviced. 

A departure of an item from one subsystem consti- 
tutes arrival of an item to next subsystem. Tho arrival and 
departure of the items are described as follows : 

Tho work station can pick up an item from 
as soon as is free to servo another item after serving 
last item it had. This is subject to tho condition that 
the work station S. does not fail in the mean time. If there 

l 

is no it on in 3. then S . waits till an item is available 

to B.. from S. and this process continues. This is valid 
3.-1 3.-1 

for all 3^, i = 2,3,,**,n. This is not valid for 
because there is always supply of items. It is ensured that 
only a single unit will be served at a time by the work 
station. 

The work station S. can. deliver an item to Bj 

i 1 

a.s soon as it completes tho service on tho item* If B^ is 
full, tho work station Si waits till an item is removed 
from B* by work station S. _ and then delivers the item 

JL L 
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to B^% This is valid for all wort stations except the 
last one, for which unlimited buffer storage is available. 

3.4 PROBLEM STATEMENT 

The problem of line balancing may be stated as 

follows : 

Given, 

1. the total number of tasks to be performed in the 
assembly line (i.e., the total work load of the line) 

2. the precedence relationships among the tasks, 

3. the mean and the standard deviation of each task time, 
and 

4. the total number of work stations available to the 
assembly line. 

Then, the problem is to assign the tasks to the work stations 
of the assembly line, following the precedence restrictions, 
in order to optimize the objective function (the various 
objective functions are discussed in the following section). 
This is to be done accounting for the following features of 
the assembly line, 

1. Variability of the task time, 

2. Failures of the work station, 

3. A case where cycle time is less than at least 
one of the task times. 
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The problem requires following data as the input to the 
model : 

1. Total number of task in the assembly line 
2* Mean and standard deviation of each task time 

3, The precedence relationship of the tasks 

4. Total number of operators available 

5* Failure characteristics of the work stations 
6, A starting value for the cycle time. 

3.5 THE OBJECTIVE FUNCTIONS 

The various objective functions related with the 
cost of the assembly line systems are : 

1* Minimization of cost per unit, 

2, Maximization of profit per unit, 

3. Maximization of profit rate. 

Although the management may be interested in the 
objective functions related with the profit functions of the 
system, yet it cannot be justified here because of its 
subjective nature. A distinct estimation of contribution 
of various productive and nonproductive assembly operations 
to the total profit earned Is very difficult. Therefore, the 
present work restricts itself to the 1st objective function, 
namely, minimization of cost per unit produced from the line. 

The following chapter is devoted to the methodology 
to solve the problem stated in section 3.4. The model deve- 
loped in chapter 4 has been validated in Chapter 5. 


CHAPTER. IV 


DEVELOPMENT OF MODEL AND SOLUTION PROCEDURE 


The prohlera stated in section 3.3 has been 
divided into two phases i 

Phase I : Balance the line -using mean task times and 

determine the maximum output rate that can be 
achieved from the line. In this phase., the 
task times are treated as deterministic. 

Phase II: Determine the optimal levels of buffer invento- 
ries between stations accounting for probabilistic 
nature of the task times and failure characteris- 
tics of all the work stations, which minimizes 
the total cost of production per unit produced 
from the line. This phase also determines the 
expected output rate from the line. 

Both the phases have been discussed separately, 
in the following sections. 

4.1 » 

4,1*1 Assumptions s 

The following assumptions are made to develop 
the model for this phase : 

1. Only one operator is assigned to each station, 

2. Ho operator is fixed for any particular job# 
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3. T-To operator is absent. This nea.no that the management 
ensures that the needed number of operators are 
available on the working days. 

4. It is a single model assenblj* lino. 

5. There are no zoning constraints. 

G. The attrition allowance, i.c . , the percentage defective 
is independent of the output rate , 

7. The task times are independent of one another. 

8. The efficiency of all tho work stations is same. 

4.1,2 il omo nc 1 atr.ro : 


1. 

m 

• 

• 

total number of tasks in tho assembly 

2. 

n 

* 

• 

total number of operators available 

3. 

ih 

» 

• 

mean task time for i"*-'* 1 task, i = l,2,.,,,m 

4. 


• 

« 

standard deviation of task time. 

5. 

H 

O 

« 

* 

initial cycle time 

6, 

F. . 

« 

• 

3"^ element of i^* 1 row of Immediate Follower 




Matrix F. . indicates the task number, which 

13 




immediately follows i^ 1 task. 

7. 

NT., : 


total number of tasks immediately following i 




task 

8 . 

n* 

* 

* 

maximum of ITT., i = l,2,...,m 

9. 

G. . 

1 m 

• 

• 

3’"^ element of i^ row of a binary r-iatrix 


called f Follower Matrix 1 

G. . = 1 implies that 3^ task immediately 
J 

follows i^ task. 
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Cr n - . = 0 implies that 3 ^h task does not 

J 

immediately follow i "^ 1 task 

10 , H.j s 2 ^ olomont of i ^ 1 row of a "binary matrix 

J. J 

called 1 Total Follower Matrix r 

Hj. . =1 implies that 3 task follows i^ task 

J- J 

IL j = 0 implies that 3 task c"oes not follow 

— J 

i ^* 1 task 

11, t Cycle time for k "^ 1 iteration 

“f“Th 

12, If. i i element of the sequence of tasks generated 

"by the method of line balancing 

13, SI : Slack time left at work station^ 3 =l,,,,jn 

3 

14, t • : number of tasks assigned to 3 "^ station 

J 

15, djr : balance delay for k^ iteration 

16, T. : Time available with station till latest 

J 

assignment of task there, 

17, A_- , B, and D, : are arrays as discussed in section 4,1,4 • 

1 5 1 1 J 

IS, th , ik g , ¥^ 3 , : are weight ages as discussed in 

section 4,1.4 

IS, Rj : Summation of task times of the tasks following 
i"kh task. 

Few other notations for the chase are described 
later in this chapter as the model is developed, 

4,1,3 Balancing the Line : 

The implication of balancing a line is to assign 
the various tasks of the assembly line among the various 
work stations uniformly and compactly while following 
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the precedence restrictions. Therefore, the best balance 
is one 'which simultaneously achieves maximum compactness 
and maximum uniformity in the allocation of total work load 
among the work stations. Balance Delay Ratio (BDR), a 
measure of total idle time of the line, is one of the most 
commonly used measures to represent the compactness 
aspect of work load allocation among the work stations. 

It is defined as follows i 


m 

n x - C M. 

i=l 1 1 


BDR = - 


n x c- 


where n, c,^, P'1 are as explained in section 4,1*2. 

The smaller the value of BDR, the more compact 
is the assignment, A zero BDR value implies that there 
is no idle time in the line and the line is perfectly 
balanced. But for most of the practical problems the 
minimum achievable BDR is a non-zero quantity because of 
the descrete values of the task times. This indicates 
that some work stations -are more heavily loaded than the 
others and that the line contains some idle time. It is 
desired then, to distribute the total idle time evenly 
among the work stations and make the work load allotment 
more uniform* 

Almost all the line balancing techniques, repor- 
ted in literature achieve the first objective only i.e«, 
to attain the minimum idle time of the line. Uniform 
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distribution of the total idle time among the work sta- 
tions is not been reported so far. This objective can 
be achieved if several solutions having the same minimum 
value of ' he BDR for the balanced line are available. 
These solutions are then used to select the one which 
gives the most uniform distribution of the total idle 

time among the work stations. 

The problem handled In Phase I Is divided into 

two parts to achieve the following two objectives s 


Objective 1 


Objective 2 


To achieve maximum compactness of the alloca- 
tion of tasks to the various work stations, 
that is, to achieve minimum BDR. 

Having achieved the balance with minimum 
BDR, to distribute the total idle time 
uniformly among the work stations. 


4.1.4 Part I i 

The aim in this part is to accomplish objective 1 
i.e. to achieve minimum BDR. COMSOAL technique develo- 
ped by Arcus [l] has been chosen to achieve this 
objective. The reasons for this choice are given below. 

1. Master and Others [27, 17) have reported that COMSOAL 
technique is found to be the best for large size 
problems (problems having more than 40 tasks). 

They claim that this requires much less computer 
time as compared to other techniques e.g., the 
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Dynamic Programming approach by Held (13). It also 
yields optimum solution in most of the cases, 

2, This technique, unlike other techniques gives 

more than one solution which have the same value of 
the BDR* These solutions are then used to achieve 
the objective 2, as will be discussed later. 

The essence of COMSOAL lies in generating the 
feasible sequences of tasks (i.e. the sequences of the 
tasks obeying precedence relationships) and to assign 
them to the various work stations. The sequences giving 
minimum number of operators are said to form the list 
of best ones* To achieve the maximum output rate for a 
fixed number of operators the same procedure is repeated 
by reducing the cycle time gradually till the required 
number of operators for a cycle time exceeds the number 
of operators available. Thus, the sequence giving 
minimum cycle time for a gj.ven line length are the best 
one s . 

Since the generation of every feasible sequence 
for a large scale problem requires enormous effort 
Arcus suggested a few heuristic approaches to generate 
the feasible sequences. This approach reduces the compu- 
tational effort and at the same time gives optimal solutions 

the 

in most of the cases. Among/many heuristics suggested 
the following have been selected in this work which were 
found to be the best by Arcus, 
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1. 


Assign weightagcs to the tasks in proportion to 
their moan task tines* The weightagcs are calculated 
as follows : The weight a go for i ul candidate task, 


\1 


il 


Ml 

1=1 


where, is the total number of tasks which ate candida- 
tes for assignment to the work station being considered 
p re sently* 

2. The tasks are weighted in proportion to a number 1/X 
whore, X is equal to the total number of unassigned 
tasks less the number of all the tasks that follow 
the tasks under consideration less unity. Fathematica- 
lly, the weight ages can be represented as follows : 


b, e 
1 


• 1 


(k-fi-1) 


w i2 = X 


i=l 


b i 


whore , 

f^ s total number of tasks following i bb task 
k : total number of unassigned tasks at the time when 
i bb task is candidate for assignment. 

3. A task is weighted in proportion to a number deter- 
mined by adding 1 to the total number of tasks 
following it. Kathematically, it can be expressed 
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as , 


W 


13 


(£ + l) 

~Y 

Z(f i+ 1) 

« “t 


4: » 


A task Is we I ghtod in i^roportion to the 
of task times following It and including 
The exnression for the weightago for i^ 


task Is, 


W I4 = 


( % + _ R i) 

Zc v 

i=l 


summation 

itself. 

candidate 


It is to- bo noted that all the weights are normalized. 
After all the above mentioned weightages are determined, 
final weightages W n - are calculated for all the candidate 
tasks by multiplying different weightages of the same 
tasks. That is, 

\I. = X VI j. g )<. for, I = 1,2,..,,^ 

These final weightages W. are also normalized to 
obtain an array r n - of such final normalized weightages, 
for all the candidate tasks ready for the assignment* 

4, 1.4,1 Generation of feasible sequences for balancing 
the line - 

Following are the steps needed in generation 
of the feasible sequences of the tasks. 
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Step 1 : 


A matrix F. . is formed such that the elements 

— J 

til 

in i row of this matrix indicate the tasks 
for which task precedes immediately. 


Step 2 : With the help of the matrix F . . . form an. 

a 3 

array ■ for all the tasks of the assembly 
line such that i^ element of this array 



represents the total number of tasks immediately 
preceding i u task. 

Form a set containing the tasks from array 
A i which have zero value in array Aj_. 

Scan the tasks in the set satisfying the 
following condition : 



Here 3 is the station being loaded at the time 
of latest assignment of task. Place all such 

tasks in the set D . . 

1 

If no such task is found in set B. for 
which above condition is satisfied, consider a 
next work station with the new value of 
Tj = for the fresh assignments of candi- 
date tasks in set D. * Here c-^. Is the cycle 

-j* Vi 

time being considered for the current k 
iteration. Step 4 ' is then repeated with 


the new values 
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Stop 5 : All the elements of set D. are the candidate 

' " f '*“ 3, 

tasks for the assignment to the work station 

under consideration. All the tasks in the 

set are now weighted according to the 

procedure described in section 4.1,4, The 

selection of task from set is biased in 

accordance with the weightages assigned to 

the tasks and is not purely random. The 
th 

g ts.sk selected in this manner is allocated 

4-V, 

to the j work station. This allocation 

results in the reduction of the time available 

to this station by an amount jl • That is, 

th ^ 

time available at i station for other 
allocations is now (T^-/^). The number of 
tasks t. assigned to station j .is also 
incremented by unity. 

Step 6 : After the q^ 1 task is assigned to station j in 

Step 5 the values of A^ are modified as 
follows : 

A = a high value, say 10,000 

hi 

i Vj e 

Aj = A-jkj e F g< ^ 3 * 9 
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Step 7 _ : The steps 1 to 6 are repeated to generate a 
sequence of tasks till all the tasks are 
assigned to the work stations. 


A number of such sequences are generated by 
repeating the above procedure. The sequences- 
yielding minimum number of work stations are 
selected. 


The urocedure described so far does not include 


the case whenjkjL c^ for at least one i. The following 
section describes the procedure to account for such a case. 


4. 1.4. 2 Consideration of a case wherejJL. c^. for 
at least one i . 

As discussed in section 3.2 a general model 
of line balancing should indlude this feature , To take 
this case into account, the procedure a.dopted in this 
work is to introduce parallel stations. Since the task 
r for which |JL^> c, f . cannot be further subdivided into 
smaller tasks, the station time of each of the parallel 
stations must be greater than the value |<l r . In fact, 
all the parallel stations according to the procedure 
adopted here will have work station times as the integer 
multiples of the cycle time such that the task r can be 
performed within this higher value of work -station times 
for the parallel stations. All the work -stations are 



32 


considered exactly identical as far as the nature of 
work load assigned to them is considered. If all the 
para3.1el stations are taken together and are considered 
as a single equivalent work station the effective cycle 
time for this unit will be equal to cycle time for the 
rest of the stations. 

After assigning task r to the parallel stations, if 
any time is left with them, it may be used to allocate more 
tasks to the stations, if the precedence restrictions permit 
to do so. The details of the logic used for this case are 
given below. 

The number of parallel stations required to perform 
the -task r 


P r = 



where, f is a fraction (0 <Cf<'l) added to make p r an 
integer. 

Thus, effective time available with all the parallel 
stations = p r x c^. After assigning task r to parallel 
stations, the time available with each of the parallel station 
fob assignment of other tasks is. 


P* = P„ K 
a r r 
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For the task r a new term * Virtual Task Time 1 is defined 
here as follows : 

. t 

Virtual Task Time U- = c,- - p f 

' r r 

r 

It Is to he noted that U is not the actual task time for 

r 

task r. The balancing is then carried out according to 
the procedure discussed in previous section, with the 

t 

difference thatji is now replaced by jS, and n is replaced 
by a new value of n = (n-p, +1)* When the final balance 
is achieved, the station having ta.sk r is replaced by p r 
number of similar parallel stations, loaded identically* 
The flow chart in figure 2 shows in logical 
form the procedure developed for attaining the first 
objective. The nest section explains the procedure for 
accomplishing the second objective. 


4.1.5 Part II : 

In this part the task sequences obtained in 
part I are examined again with a view to uniformly distri- 
bute the total id? o time of the line among the work stations. 

To evaluate the uniformity in distribution of the 
total idle time among the various work stations, an index 
called the "Smoothness Index" has been introduced. Smooth- 
ness Index is defined as, 

/ m 


= Y~ (s< 

/ i=l aV 


- St )‘ 


S.I 
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C Start ) 
=jF— 

f Head the Input "\ 
\ d ata j 


Calculate theoretically mini- 
|mum possible cycle time' c m j n 

m 
c 


mm 


IfV n 

i=l 


lake initial cycle 
time 

1 r 



Optimal solu- 
tion has 
reached 



Yes ^^cyclo time takeri\ 



Find out the 
number of para- 
llel stations 
required 


mo 


Find out the number 
of ope r ato r s , now 
available for rest 
of the assembly 


T 


Ilodify moan task times 
for these tasks 


Balance the line 
using COMSOAL 9 as 


discussed in 
sec. 4.1.4 


n is 
smaller 

<c 



Compare .. \ n is 

minimum number of \ grea- 
operators required Vie r 
by above balance 


with n 



© 



both are 
equal 


(CONTD...) 
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where, S.I. = Smoothness Index 

S« = Slack time at ith work station 

S av = tlie avera S e of the slack times of all the 
work stations 

S.I. is a measure of the distribution of the total idle 
time of the lino among the work stations. The interpreta- 
tion is that for the task sequence having the minimum 
value of S.I. , the dispersion of the slack times of the 
work station about its mean value is minimum. This suggests 
that the slack times are distributed- most uniformly among 
the work stations. In fact, the smaller the S.I. more 
uniform is the distribution of the idle time. 

The task sequences having S.I. values close to 
the minimum value of S.I. may be treated as equally good 
sequences, and thus provide the line supervisor with 
alternatives. The line supervisor can then select one of 
these alternatives according to his own convenience. 

It should be noted that the definition of S.I. as 
given above is different from the one proposed by Hoodie 
and Young |28j • They have defined S.I. as, 



where, S. - station time of i^* 1 work station 

7 l 

S 

max 


maximum station time in the line 
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They have used the index in a different context* They 
have employed it for trading and transferring tasks from 
one work station to another for the purpose of making work 
station time nearly equal. They have not considered the 
compactness aspect of the allocation of the tasks to the 
various work stations. 

The procedure of balancing the line, as explained ' 
in Phase I has been computerized. Every care has been 
taken so that it can be applied to any general line balancing 
problem. The listing of the computer program is given in 
Appendix B, The comment cards have been put liberally to 
explain the working of the program and meaning of various 
important variables used in the program. 

Phase I is followed by Phase II , which is explained 
in the next section. 

4.2 P HASE II : 

Phase I yielded a balanced lino, where the work 
stations are loaded with tasks uniformly and compactly. 

This was done using the mean value s of the task times. 
However, In real life situations one comes across 

1. Variability In the task times ,and 

2, Failure of the work stations. 

The objective in this phase Is to account for 
the above mentioned two aspects by incorporating optimum 
■buffers between the work stations. The values for buffer 
capacities are determined with a view to minimize the total 



cost of production por unit produced from the lino. 

An assembly line represented schematically in 
figure 1 can be considered as a single channel queuing 
system with n subsystems each consisting of a service 
facility (work station) and the item (the buffer inventory) 
in the queue waiting for the service. The service time 
and service faility 1 s failure are random phenomenon and 
have some probability distributions associated with them* 

The assumptions listed below have been made to 
develop the model. The basis of these assumptions is the 
characteristics of the various types of data collected for 
the assembly line of J.K. Electronics. The details about 
the data collection will bo discussed in Chapter 5, Minor 
'modifications in the methodology can incorporate different 
probabilistic characteristics of the data. The assumptions 
made arc as follows : 

1. The time between two failures of a work station 

exponent ially 

is a negative £ distributed random variable. 

2. The repair time for each station is a normally 
distributed random variable, 

3. The service time of the work station is a 
normally distributed stochastic variable. 

The mean' and the standard deviation of the service 
time of a work station are given by the following rela- 
tionship s : 
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whe re , 


H . 

S 1 


h 


1 


mean of service time for I^ work- station 
(also called work-station time) 
mean of the task time for task 

set of tasks which are assigned to 
work station in Phase I 

r l,^,..*,n 



whe re , 

CT : standard deviation of i ^ 1 work station time 
standard deviation of the task time for 3 ^ 

t3 

task 


It is to he noted here that the task times are 
normally distributed variables. Also to ho noted is that 
in the expression for standard deviation of the work 
station time, the correlation among the task times of set 
k% has been neglected. 

The transactions in the assembly line being consi- 
dered were explained In section 3.3. This assembly line 
with the features listed about namely the variability of 
the task times and the failure of the work stations cannot 
be solved analytically because of the complexities involved. 
Therefore, digital simulation has been resorted to solve 
the problem. A logical model Is developed in the present 
work to represent the above mentioned assembly line system 
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1, n 

k. h* 


v, 


S3. 


’ * Mri 


to to solved on a digital comnu'cer. 

4 .2.1 he rminolog-’ : 

In-'iit data to tlio system : 

total number of work- stations in the assembljr 
mean work-station time for i r,h work station 
standard deviation of i^ 1 work- station tine 
parcneter of negative exponent tally distributed 

-}-i n 

intorfailure tine of i l,1 ~ work station 
5 : near, and standard deviation of repair tine 

for i 1 '* 1 work- station 

is the number of stations in the group of 
is the member station of j group of 
is the number of items produced after which output 
rate is computed 
Other important variables t 

S. aj : indicates the value of time of the next arrival in 
work station 

10. : indicates the value of time of the __ next departure 

ioli 

from the i work station 

11. f . : indicates the value of time at which i work 

station fails 


6. g i 

7. u 


•3 


8 . 


nr 


1 


12, q ; indicates the value equal to the length of queue 
after work station ( i . e . 3 number of items in 
buffer i) 
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indicates the status of i^* 1 work station 

ST. = 0 imnlies that the station j. does not 
1 

have any item under process and it is 
ready to tale one item from the buffer. 
•3T.=1 implies that station i is now occupied 
hy an item. It irj.ll deliver the item 
as soon as the service is complete, 
has the value of time of the most recent change 
in the status of the stations. It is some sort 
of clock time which gives an indication of the 
time the process has already taken place, 

4-Vj 

denotes the number of items produced by i 
work station till T, 

4"Vi 

denotes the units of tine for which i buffer 

space has (j-l) items in it 

for i = l,2,...,(n-l) and 

j = 1,2,... , (Bj.+ 1). Here is the 

th 

buffer capacit3 r for i buffer space, 
actual time taken by i^- 1 work station to produce 
an item 

has the value of the time of the previous arrival 
in the i"^ work station. 

has the value of the tine of the previous depar- 
ture from the i i:h work station. 
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Following arc the design variables of the 
system : * 

^ : contains the value of buffer capacity 

for 1 th group of Y ± as explained in 
section 4.2.2. 

21. B^ : buffer capacity between i^ and (i+l)^* 1 
work stations. 

The other variables and parameters related with 
cost model have been described in section 4.2.6. 

4,2.2 kethod for reduction of combinatorial size of the 
Problem ; 

As explained in section 3.2, the consideration 
of different buffer capacity sizes at different work 
stations makes the problem highly combinatorial and there- 
fore practically unsol vable. A methodology has been 
developed, in this research work to unriddle this problem. 
The methodology makes possible the allocation of different 
buffer levels at different buffer spaces. 

Consider the assembly line shown in figure 1 of 
chapter 3. The buffer capacity is dependent upon the 
variances of all the stations preceding and following 
it. However, a buffer is mostly effected by the stations 
adjacent to it on either side. Further, any contributions 
due to other stations is small and therefore can he 
neglected. The magnitude of their effect is further 
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reduced duo to buffer allotted between thorn. Therefore, 
the buffer capacity is mainly dependent upon the variances 
of stations and S_. + ^, Following this assumption , 
associated with each is a number Y. to help in the 4ev&lopmen- 
of the model. The number V- is defined as follows : 


■yn = Variance of station + variance of station S n - + i 
The following scheme is followed to reduce the 
combinatorial size of the problem. The steps are given 


below. 


1. Arrange all the i = in the increasing 

order of the values of V. 

i 

2. Divide the whole range of V. into k ( k <C n-1 ) subdivisions 
such that each subdivision have equal number of buffer' 
spaces. If it is not possible to have equal number of 
buffer spaces in each of k groups then the buffer 

say Bjj is included in the group to which V^ is 
nearest. 

3. Assign sane initial values of 3.. to all i belonging 

to group 3 . Call these values B* , 3 = l,2,..,,k. The 
values B T . may be different for -each 3 , 3 = 1,2, 

It is also to be noted that during the simulation of 
the assembly line, whenever the level of buffer Bl 
is changed for a particular group 3 , all the 
contained in 3 ^ group are changed accordingly. 
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4.2.3 Optimisation Procedure : 

Aj.ter following the scheme explained in previous 
secu.ion, there are k (k n-1) design variables, namely, 
the values of B* for k groups of V, . The values of B* 

a - a i 

range from zero to infinity and therefore it is not feasible 

to calculate objective function for each possible combina- 
tion of their values. In order to determine optimum values 
of , 1 = l,2,...k, and also to reduce search Univariate 
Search technique [s] is used. The reason for this choice 
is r.'. inly to have economy in computer time. Most of the 
optimization techniques, other than Univariate method 
usually require local exploration at intermediate stages of 
process.’ o.g., evaluation of the gradient of the function 
which need many function evaluations. In the present 
problem each function evaluation requires one simulation 
run of the process and hence a large computer time is 
consumed. The Univariate technique on the other hand, 
does not need, this information and therefore, less computer 
time is required. Another reason for this choice is its 
simplicity and ability to give good solution for a problem 
with small number of variables. 

Fig, 3 shows tho procedure for univariate search 
schematically with tho help of a flow diagram. Following 
are the steps involved in this procedure. 












Lie optimal solution has 


reached. Last design 
point was optimal value 
of the design vector 
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1. Choose a starting vector B^ , q = l,2,..,,lr 
° • Dsternj.no wnethor the predetermined step length t should 
he positive or negative* In other words, for a particu- 
lar direction of movement, S , does the function 

4 

oe crease in tne posi.tj.ve or negatj.ve direction of s q ? 

B. Tahe a step of + t in the direction S n (whichever 

step length minimizes the function value ) and calculate 
the function value F t at the new design point B^ +1 = 

-q ±tSq. If the value of Fq at the new point is 

loss than the value F q , increase the step length t ty 
a factor e, where c^> 1, and repeat step 3 and 4, 

Qt; erwj.se set the design point B* = B ! + t„ S , where 

q+l q — p q 7 

tp is the last 1 successful value 1 of the step length. 
Choose the next S q and repeat steps 2, 3 and 4 for 
the new direction. 


4. Stop vheo no direction S , produces an improvement in 
tire •■'‘unction. 

The directions of movement, S q f s, are chosen as a 
cyclic ordering of the unit vectors e.g M Sq = (1,0,0,.,, ,0) , 

Sg — (0,1,0, » » , , 0) ,«.« , Sp- — (0,0, , ,0,1) , (1,0,0,.. ,0) # 

It is to he noted that this procedure gives only 
the local minima of the function. Therefore, it is 
preferred to evaluate local mininas from different star- 
ting points and. select the minimum of all these minlmas. 
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4.2.4 Simulation of the assembly line : 

A logical model of an assembly line which rep re** 
sents both of the following features, namely, 

1. The variability of the task times, and 

2. The failure of the work stations, 

has been developed in this research work. Fig. 4, 4(A) 
and 4(B) show the logical flow diagram for the proposed 
methodology. 


4. 2, 5. The Steady State of the Assembly line : 

During the simulation run the assembly line is 
assumed to have reached the steady state when the rate of 
production from the line is stablized. The rate of 
production is computed in consecutive periods of time in 
which t"o .line produces m f additional products. That is, 
fne ra.te of production is computed when the line has 
produced the (m*)^, (2m 1 ) th , (3m l ) th .... product etc. 

The value of m f is selected arbitrarily and j.s kept close 
to the average number of products produced in one shift. 
The rate of production per unit time is defined as follows: 


0 


m 1 


i+1 


Time required to produce (i+1) x n 1 
products - time required to produce 
i x m r products 


where, 0 ±+1 


m* 


rate of production in i" 1 "* 1 period 
a number as discussed above. 




Y arrival time is 
j minimum 



* These boxes are explained in Fig. 4(A) and 4(B). 







Figu re 4(A) 
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C Start ") 



** The stations can nick up an unit immediately after 
it has finished working on one. Therefore, next 
arrival times will he equal to current departure 
times. 


Figu re. A(BJL 
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Since the transactions in the assembly line 
are probabilistic in nature, the chances are that the 
computer 1 values of 0 3 and 0 i+1 are close to each other 
even when ino l.ine has not reached the steady state. 
This is depicted clearly in figure 5. 

The following procedure has been adopted in 


this wor ]: to avoid above mentioned situation, 

1. Let the process continue till p periods. Here p is 
the number of periods after which the assembly line 
is o::uocted to reach the steady state. Thus number 


is estimated by tailing a trial simulation run of the 
system. 

2. Calculate m 0 -pcriods moving averages of the rate of 
production in the p^ 1 and (p-1)^ 1 period, as given 
below. 


m.. 


I = i. 

p m 


o 


H o-- 

3=1 


P~i 


, _ 1 

A = — 

P-J- m. 


m o " 1 

1=0 


0 . 

P-3 


where, A^ = moving average in p^ 1 period 

A moving average in (p-1)^ period 

m Q = an arbitrarily chosen number such that 
1 <m Q < p. 
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O 

o * 


Go v :* ■ * 1^0 xi and * I f 
than (a s-all quantity 
desired) tho steady?* state 


their difference is less 
5 depending upon the accuracy 
of the process has been 


reached. Otherwise increase p by 1, simulate the 
process for the rent period and. repeat steu 2. 


taking m 0 ~periods moving average 


smooths out the 


randomness of the values of 0^ and therefore gives a 
better indication of the steady state. 


4.2.6, Stochastic Convergence of the Rate of Production 
from tho Line : 

A tost for stochastic convergence is conducted 
after tho assembly line system reaches the steady state. 
For this mirposo the sample size is increased by increa- 
sing the total length of simulation run. The sample size 
is calculated for the desired level of significance. The 
details of the procedure are given in Chapter 5, 

4.2.7, Variance Reduction in the he a sure mo nt of the 
Production Rate : 

Following measures have been taken to reduce the 
random errors in the measurement of the production rate, 

1, The use of sane set of random numbers for 
different simulation runs leads to reduced 
random error for comparison of the stochastic 
variates. The values of the stochastic 
variates obtained by using same set of random 
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numbers are lively to be positively correlated 
and therefore the me a sure me nt of the difference 
of the stochastic variates will have reduced, 
random errors, 

based upon the above principle. the same 
sot o" rend om numbers have been used for 
si: aviation runs for different design points ■ 
in the optimization procedure. 


2, The precision of a given sample is improved by 
using the concept of Antithetic Rule [34] . 
According to this rule, two simulation runs, 
instead of just one which is done usually, are 
taken. The simulation runs are identical, 
apart from the fact that the random number 


sequences used therein are different in both 
the runs. If x_ , x 2 ,...,,x n is the sequence 
of random numbers used in the first run, the 
set to bo used in the second run will be 


(l-x 1 > , (l-x 2 ) , • ••( 1-x ) • The justification 
is that these two sets of random numbers are 
strongly negatively correlated and therefore, 
tend to produce results on the opposite sides 
of the population mean, and when taken 
together, will give results closer to the 
population mean than would be likely o the wise. 



A 0 


.8. The On «;+. bo del : 
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...io costs have ‘been Included into the 

total cost of ■-•roc net. Ion : 


3.. Fixed labour cost 
8. Incentive labour cost 
£. buffer lidding cost 
•1. Fixed overhead cost 
•1 . Variable overhead cost 
S . ' "ato ri al cost 


1. Fixed Labour Cost : 

Total fixed labour cost/ shift, L c = n x x^ 

■where, n : total number of opera- 
tors in the line 
x n : labour wage per labour 
per shift, 

2. In.con.tivo Labour Cost : 

Total Incentive cost uer shift, I = (0,, - CL.) x x Q x n if 

* Q b III d, 

°s> °m 

= 0 otherwise. 

Here, 0 : steady state output rato per shift 

O m : minimum output rate, over which the 
incentive is paid 

x : incentive cost of producing an 
additional item. 
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Bu f .'V: t ■ ho? <1 Cost : 


n 


Total 


i'.T : ‘or holding cost per shift, H c = } q 


i=2 


x x- 


hi 


’ *ho ro . o 


, Q. : average buffer inventory level at 


‘hi 


Fired Oyn a- ad 


. th . „ „ 

i Duffer space 

Cost of holding buffer at 1 th buffer 
space/ shift 
Cost nor shj.ft = F. 


5. Valuable Cvo head Cost : 

Total variable overhead cost per shift, OV - J z (I +L ) 

c c 

-'her-' ,7 is a specified factor relating overhead 
or.' si*, vf th total, labour cost. 

6 . I late : -j. n.l C <“■ : ; i" : 

Total r a to rial oort per shift, li c = 0 g X x n 

where, is the rate rial cost per unit 
The total cost of crodurtion por unit produced 
from the lino in obtained by summing the six costs listed 
above. Fat horn t; i c ally, the total cost of production per 
unit .is err rosso cl as, 

OP = ~~ I I + I + H + F c + OV + ;i 

°s L c c c 

The other important parameters to be computed 
during simulation are given on next page. 
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1* Force nuayo j 1 3.13. z at .ion of the work station. 

Force iv, utilization of tlio work station means 
the nor cent of t'-e total time the operator is busy. It 
is defined as follows : 

IT 1 

PU, = •--- 

J. n 
JL 

’-hero Plh : percent utilization of the work 
station 

TJj^ * total productive time of the work 
station# * 

I : total clock tine of simulation run# 


2. Average buffer inventory level at i^ 1 buffer space* 

This .is analogous to the average queue length in 
a queuing system and .is defined as follows : 


i 


if 1 x 


T 


where : average buffer .inventory level at 
fh 

i buffer space 

f To 

Fu : buffer capacity of I buffer space 
Q 3 - as defined in section 4.2.1 
The methodology dove-3 oped j.n this phase has also 
been computerized. The listing of the program Is griven in 
Appendix B. By making minor changes in the computer program 
one can Incorporate probability distribution for the inter- 
failure and repair times, other than the ones taken in 
this work. Tho next chapter discusses the data acquisition 
and the validation of the model. 



C JayPTER V 

DATA ACQUISITION .HD ANALYSIS 

The model developed in Glia;) tor IV has been 
tested on an assembly line of I. II. Electronics, Kanpur, 
one oi' the load in a TV set r nnufactu.ro rs in India., 

T'.'o assembly line studied is called the 1 VIF-SIF 1 
assembly lino. This is the main subassembly line of J.K. 
Electronics. The output from this line, VIF-SIF plate, 
is one of fir? major components of a television set. The 
line eon, r lots of sin work stations and each work station 
is manned by an operator. The management wants to retain 
the nurd or of operators as six and wants to improve 
tho product:' v\ty of the line by using scientific planning 
methods . 

Hi the VIF-SIF assembly line each operator is 
supposed to perform a prescribed sot of tasks. Due to 
the pro sene? of human element, the time contents of each 
task are probabilistic in nature. Therefore, the time 
needed for a set of "asks to bo performed on a work station, 
is also probabilistic in nature, A time study was 
conducted to dotormi.no the exact nature of the distribu- 
tion for tho tas v times. As will be seen later in this 
chapter, tho task times are- normally distributed. Using 
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tine there ’Jill l.-o no output from tho work station, this 
period has ’con eons id?’ od as tho failure period of the 
work station. The other failures o . g . , power failure 
etc. are not, considered. This is because these failures, 
sort of, freeze the process of assembly lino altogether. 


They stop the working of the whole of the assembly line 
anc 1 tho '’c fore, neither they affect the general nature of 
the flow of tho items in the line, nor they affect the 
different, levels of the buffers. 


The phenomenon of failure of work stations and 


the repair time (i.e., the time duration after which the 


operator comes back to his position and starts working 
again) arc also found to bo probabilistic in nature. As 
will be seen later in this chapter, the interfailure time 
of a work station follows a negative exponential distribution 


Tho Central Limit Theorem states that if xp are continuous 
and independent random variables with mean value r-p and 
standard deviation value , then tho sum x = xi+xgH- • • «+x^ 
will approach a normal distribution with mean, value 
p> = j^' + M 2 + , . , , + fife and standard deviation value 


O' n Sj g 2 + 


+ O' 


2 

k 


* 


* 9 p * 
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and tr. "• repair time 
normal di str.i butinn 
i.n the S'lct.io.: 5.6, 


of the work st.at.ion follows a 

The details of this nj.ll be discussed 


set 1.031 describes the total work load 


t’ 10 


1 T? 


Ubl 


.1.7..-1C , 


5.i r.n vci.il lc.jd i: t to-sif assefely life 

Th-'* current procedure adopted for assembling 

tho VIF-S1F -'date involves tie following stops for each 

of two operators : 

1. Pic]: up o. fresh plate fron the buffer, 

2. Fix a known nurber of the components (as specified by 
the management) onto tho plate , 

3. Flip tho plate, 

4. Bring tho soldo ring rod to tho plate , 

5. Solder all tho terminals of the components fixed in 
stop 2, 

6. Put the soldo ring rod back to its place, 

7. Bring tho pliers to the plate, 

8. Brim tho terminals of tho components soldered in 
stop 5, 

9. Put tho pliers back to its place, 

10. Flip the plate again while placing it in the next 
buffer after completion of the 30b. 
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In t: o o.l') o vo lint stops 1, 3, 4, 6, 7, 9 and 10 
fom the t a sirs which arc exactly identical for all the 
work s t o is and arc performed once in each cycle • The 
tael's in store 2, 5 and 8 arc different for different 
nor: static. >c arc! aro discussed below. 

dtop 2 cc rip rises of several tasks where each 
tasl: cons-’ sts of fixing a component, or a set of components, 
onto tk? elate, The tasks of fixing the components onto 
the plate fellow a Wnown precedence restrictions ho cause 
of t to technological constraints. Fig. 6 depicts those 
precedence restrictions. In this figure the i^ 1 node of 
the precedence diagram represents the fixing of the component 


The significance of the numbers shown on the north-east 
corners of the circles representing tasks will be described 
later in this chanter. 

The task of soldering all the terminals of the 
components f.ixod at a work station (refer step 5) comprises 
of several tasks where each task consists of soldering the 
terminals of a component. The following relationship is 


assumed tc hold good. 


n f 

°12* ...n 


si + si, + 


* * • • 



where, Si 0 .. total soldering time for soldering the 

terminals of components l,2,...,n, 
all cone together. 

s! : soldering time for soldering the terminals 
i 

of i'*'* 1 component only. 
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The above 
•it i s arc c.dr .i.ti ve 


relationship suggests that the soldering 
in nature. The sahc type of assumption 


is r ! £f’o for t rimming also (refer step 3) i.e., 



= «T»! 

. .n 


+ T l + 

2 


• * m • • 


+ T t 

n 


re 


T 1 


Tl 

1 


= total trimming time for trimming the 

.n 

to r? ■ inal s of the comp one nt s 1 , 2 , . . , , n 
all done together 

= trimming time for trimming the terminals 

4* In 

of i vi component. 


Table I shorn the cl 3. stribution of the tasks among the six 
operators os per the practice currently being followed at 
j.K. Electronics- The tasks which are exactly identical 
fo?- nil ' in work stations, are listed at the end of the 
table. Those tacks are not included in the precedence 
diagram. This is because the sum of times for the tasks 
common to nil the work stations can be directly subtracted 
from, the total cycle time of the line in order to find the • 
time available for the allocation of the other tasks to the 
work, stations. 

The following notations have been used xn the 

table : 

P : Bend the terminals of the component by 
ih 

hand before inserting it. 

P , : Bend the terminals of the component by 

oh 

hand after inserting it . 
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^Itj 1 i;Cilu too te minals of the component by 
pliers before inserting it, 

Poo : j -' on -i too to m Inals of the component by 
pliers a^bor inserting it, 

*3-^ • Straighten the component by pliers 
I ; - - 'insert the component in the specified hole 
in the nle.te. 
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TABLE I 


fcVrro s- 0 
8pondJ.n:;5 

jlo. jcorino- 0 


|numhe rs 

ITIX'-’I ”1 

tation lio.l 


Task Description 


5 Mean 5 Standard 0 

OohservedjiDeviationO 
5 task ti- 0 ( in seco-0 IK' 
One (in Onds) 0 7 

0 seconds).^ ( 

. £ 4_ 0 5 _3l_ _6_ 


0 . 


1 

'1 rite the identification 
number on the plate 

15,45 

2.31 

7. 29 

2 

Pick up 3...3 Il.He^st.ers (3)* } 
P ih’~ n 

25.47 

3.35 

8.13 

3 

P 1 ck up 4.7 Registers ( 2 ) , p . , 
I n 

13.12 

1.S7 

8.24 

4 

Pick up 680 lhRe_gj.sters(3) , 

hh> in 

19.23 

2.55 

10.77 

5 

Pic]' up 1.0, h Re piste rs (3), 
ih’ n 

13.59 

2,90 

3.20 

6 

Pick uu 5 60 ih . 2K v 5 6K 

Jlojimqrs., P. h , I n ' 

18.90 

3.17 

0.34 

7 

Pick up 10K Eejgisters„(3) , 

P ih 5 P n 

19.85 

3.19 

1.62 

- 

jjl* 

Solder all the terminals 

(36) of station 1 components 

121.51 

13.50 

0.96 

- 

T rim* * * all the t e rrninal s ( 3 6 ) , 
of station 1 components 

37.47 

4.93 

4.85 

l'Io..2 



• 


11 

Pick, up _5_G0 (2) ; 2.2 k 

Registers,' P", ? oh 

31.30 

4.19 

C C o n*fc d» * » ») 

4.59 






Ihe figure in the parenthesis indicates the 
aumher fa of similar components which have "been 
combined to form one task. 

mi t - e va] ue does not include the time for 
plciing Sp lolde ring rod and pladJag.it hack to 
its place. 

^Thi*; timp vhnp does not include the time for 
pickfng~up the pliers and placing it hack to its pla.ee* 
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TJLBLJl I__( COHTD . ) 


rp 


i 

u. 


6 


T 0 


1^ • 

14, 

15, 


17. - 

18. 

Station Ho. 3 


19. 

32 

34 

20. 

18 

21, 

29 

22, 

22 

23. 

28 

30 

24, 

- 

25 

- 

26 

46 


5 

5 


j „ 

M"*: 11 - 0.01 IT Polystor 

C1>, i -.wisut s ) , S t ,P ip ,I n VP op 

64.39 

6.16 

4,44 

Pick \’p 1_*0M Registers (2), 

s,., P. , I Vp “ 

t’ IP 3 iV OP 

22.18 

1.90 

10.67 

rick up 820 -si ,380 -d> Registers , 

3... ? p; , i;;;~p " ’ 

:4 ir n’ op 

18.50 

2.30 

10.34 

Pick up 330 -p L(2_)_3_ ,0.2 k Registers, 

P P ~ 

ih ’ oh 

18.98 

1.91 

6.18 

Pick up 68QK, 3 , SIP Registers, 
ih’ oh 

13,11 

1.63 

7.53 

Pick up 25: T Condenser. Sx., P. , 

I , p 0 IP 

n’ op 

19.07 

2.15 

8.71 

Solder all the terminals (36) of 
station 2 components 

143.82 

13,87 

2.79 

Trim all the terminals (36) of 
station 2 components 

45.67 

4.50 

1.34 

O 

r; Pick up Pin up Capacitors(2) , 

12.19 

2.49 

6.68 


* '•'M/ *yj5 

drop some molten vax over 
it ' 

Pic? i im S-Disk tj£po_ Ca pacl tpr.CU-i 

V P oh 


Pick up Dlsk_ t vp e Car ac i t o r s ( 5 ) 

T p *** 

A n’ oh 

Solder all the tenninals(36) of 
station 3 components 

Trim all the terminals (36) of 
station 3 components 


46 Pick up hack Pin jip_ c apacito r , 
I n , Solder 


17.92 3,07 0.975 

12,45 2.28 7.03 

27.51 3,52 0.48 

27.94 4.76 0.38 

119.84 15.84 1.99 

52.06 5,21 1.78 

14,44 3.00 2.38 

LSonjA*-*-*— vl. 
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TAEU3. I„(COHTD.) 


1 1 

! 2 

I 

f 3 I 

.1 _ I 

4 1 

5 

5 

f 6 

JS 

11. 

3.5 

Pi cl uio 0.01 JIF Poly ste r 

64.39 

6.16 

4.44 



Con; enscrsC 3 ) , S t ,P ip , I n , P Qp 




12, 

12 

Pick up 1.0M Registers (2), 

S,, P.“ , T *, ~P“ 
t 3 ip 3 ll 3 op 

22.18 

1.90 

10*67 

13. 

10 

Pick up G20_-^_ v 68P^^e^is.te,rs, 
b t 3 P ip 3 P n’ P op 

18.50 

2.30 

10.34 

14. 

8 

Pick up 330-^-(2j_j_ 0.2K ^egis_be_rs 3 

p p 

ih 3 oh 

18.98 

1.91 

6.18 

15, 

9 

Pick up 68 OK . 3 . 9k Re gj ste rs . 

p p 

ih 3 oh 

13,11 

1.63 

7.53 

IS. 

31 

Pick up 25jXJlQMeJl^ej^ S_t. , , 

19.07 

2.15 

3.71 



P n 3 P op 




17. 

— 

Solder all the terminals (36) of 
station 2 components 

143.82 

13.87 

2.79 

18. 

— 

Trim all the terminals (36) of 
station 2 components 

45.67 

4.50 

1.34 

Station I 

Ip, 3 




19 ; 

32 

733, Pick up Pin up Capacitors (2) , 

12.19 

2.49 

6,68 


34 

,o5 I n , Pq^ 




20. 

18 

Pick up Pe_9kinjl imQ.il_Ji ? I n > 
drop some molten uax over 
it 

17.92 

3.07 

0.975 

21. 

29 

Pick up S-Disk tynp Canaci to rills 

I , P 7 

12 .45 

2.28 

7,03 



n 3 oh 




22. 

22 

Pick up Peaking co.11 2,l n , P Qh 

27.51 

3,52 

0,48 

23. 

28 

, Pick up Disk tvne. _C a^aci_t o_r^X_5), , 

27.94 

4.76 

0.38 


30 

T P " 

n 3 oh 




24, 

- 

Solder all the terminals(36) of 
station 3 components 

119.84 

15.84 

1.99 

25 


Trim all the terminals (36) of 

52.06 

5,21 

1.78 



station 3 components 




26 

46 

Pick up hack Pin up capacitor, 

I n , Solder . 

14,44 

3.00 

2.38 






(Contd. . . . 
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TAHJIJ. (COITD.) 



Sta. 

tipn. 

Po.4 




27. 


3 Pi cl: up .Trap, coil s(3)_3_S I 
flip the plate , solder "them. 

78.06 

10.01 

1.19 


37 




23. 

24 

Pick up lOOJvPF CpndenserGU,I r ,P oll 

18.44 

2.61 

4,08 

29. 

23 

Pick up 15I-IPF ? 3 .31 CPF Co_nde_nsej?s , 
F oh 

14,33 

2.26 

6.08 

30. 

IS 

21 

, Pici: lip PP'.sJL5)., I n , P oh 

35 .95 

6*99 

1.57 

31. 

13 

Pick up Jumper 2, put wire sleeve, 

p,,, r ,T" ' 

ih ? n 5 oh 

9.06 

1.45 

4.78 

32. 


Solder the terminals of station 4 
components (ercept ITo .23) 

107,91 

15.80 

5.75 

33. 


Trim the terminals of station 4 

29. 92 

3,11 

4.32 



components (ercept lTo.28) 




34. 

38 

Pick up Preset. Potentiometer. 

S-^ 3 I, , flip the plate, solder it, 
flip rhe plate again. 

30.13 

5.45 

1.73 

35. 

17 

Pick up Jpjflpe_r _1, P- h? P Qh 

17.51 

2.42 

3.32 

Station 

HCU5 




36, 

39 

Pick up Zener. Diode, P 1h , I n ? p oll 

13.22 

1.81 

4.46 

37. 

16 

Pick up 0A7SDiode , P , I , P 

13.20 

1.90 

3.46 

33. 

26 

Pick up white capped Transistor, 
P 1h , inset it. 

37,44 

5.68 

1.72 


25 

Pick up hl-acJ:_iirAnjylsJ:nr , P- h , 
insert .it. 

15.00 

2.33 

1.08 

40 * 

27 

Pick up Integrated. Cjnicuit , bend 
its terminals by hand at a proper 

75,12 

9.12 

2.64 



angle, insert it in a particular 





direction mentioned on it. 


(Contd....) 
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TABLB I (COI'TTD.) 






41. 


Solder all the terminals (30), of 
the station 5 

91.99 

13,44 

9.34 

42. 

- 

Trim the terminals of the tiro Diodes 

8,22 

0.54 

5.64 

43. 

45 

Pick up a_2_5J!F S , , trim 

its terminals and solder it uo the 
piste 

24.06 

4.24 

6.82 

St atj.on I 

'TojhS 




44. 

41, 

42, 
43 

Pick up VIF(2), fix them on the plate 

29.89 

3.41 

8.29 

45. 

47 

Pick up SIF(2), fix them on the plate 

42.01 

4.12 

0,11 

46 » 

20 

Pick up the can and fix it on the 
plate 

19.85 

2.22 

2.62 

47. 

- 

Solder all the terminals of VIF's 
and SIF's 

163.01 

16.32 

5.76 

CO 

4 

49 

Fote down the identification 
number of the plate 

17.95 

1,96 

1. 53 

4S . 

48 

Pick up and solder the back jumper 

8.94 

1.42 

3,94 

Tasks common to all the work stations 




50. 

- 

Pick up a fresh plate from the 
buffer' 

7.12 

1,42 

4,75 

51. 


Pick up the soldering rod, bring 
it near the plate, then keep it 
back to its place . 

19.90 

4.43 

3.21 

52. 

*** 

Pick up the pliers, bring it near 
the plate, then keep it back to its 
place. 

12.19 

2.49 

6,68 

53. 

- 

Flip the plate, flip it again to 
i t s” o r i g inal p o s i t io n . 

8.57 

1.20 

9.24 
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5.1.1. 


Do termination of the tine contents for various 
tasks : 


To determine the tine contents said their varia- 
te Titles for the various tasks -performed on the VIF-SIF 
assembly line it is necessary to conduct a time study 
because tills information is not currently available from the 
factory. 

In order to have a statistically sound estimate 
of the time contents and associated variabilities for the 
various tasks, it is required that enough time study obser- 
vatiohs are made. The procedure for the determination of 
number of observations is discussed in the following 
section. 


5, 1.1,1. Statistical Determination of the Furiber of 
Readings required for the Time Study - 


m be task, times are assumed to follow a normal 
distribution. This assumption has been validated in the 
ne:;t section. The following procedure has been adopted to 
determine the number of readings for task time required for 
each task. The reader may refer the reference |23*| for 


complete details of the procedure, 

1, Choose a value for confidence level c, and confidence 
.Interval I (refer Fig, 7). The values of c and I have 
following interpretation : 

P.rob [ f - I ^ T ^ T + I ] a c 
where, T ; sample mean 

T : population mean 
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2 


m oi 


a e - nun 06 r of readings for the task. 


Calcul ate the samp le standard deviation s, as follows: 
S = 


t 2 _ U T) 2 




\! 


N-l 


Commute tie number of rocadings required N, as follows: 


~ ' V 
T !J 


XI. 


-\i~y -i-*n 

.L 


>V 


area under t-distribution curve for 
confidence level of c and degree of 
freedom of V . Here V is equal to 
value (K-l) . 


5 . If IT > It 5 take IT-*? T more readings for task trim under 
c on side ration . 

For the fine study of the VIF-SIF assembly line 
following values of parameters chosen are 

c = 0,975 5 I = 0,025 of T 5 M = 40 5 tQ.975 j3 9= 2 - 04- 

Table T 'shews the mean and th ra standard deviation 
of tlio fast", tines which are mentioned against the corres- 
ponding tasks. 

The data for interfailure and repair tines of the 
work stations were also collected. The details of the 
failure and rep-air characteristics are given in section 


5.2 
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5. 1,1. 2, Chi-Square Test for the Task Tines Distribution — 

•As mentioned previously, the task tines are 

assumed to be following normal distribution* This as sump— 

2 

tion is validated in this section using test. For 
the complete details of this test, the reader is referred 
to reference [6 ] , The computer listing, in apnendix B, 

O 

includes a general program for ^ test for normally 


c.ist ributec random numbers. The value of is listed 
in Table I for each of the tasks of the VIF-SIF assembly 

w 2 


line * 

The 

C ?.\T 

.tic a! 

97. If 

and 

cl Cl 6 

gree 

From 

Table 

i, 

one ( 

value 

s of 

n 

{ .f 

is 1< 

This 

0 o r* i 

I*T. i 3 

that 


distribution. 


5 . 2 FAILI : ITT C i ;AI ACT3RI ST I C S 

5,2,1 Interfailure time -distribution j 

Table II shows the first and second moments of 
interfailure times for each work station. In this table, 
it .is observed that for each work station, the values of 
the first and the second moments of inte rf ailure times 
are nearly equal. This gives an indication that these 
times are exponentially distributed. To confirm Chi, 
test has been applied to these variables. 


2 
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2 

Application of }jf test needs the values cf the 
parameters of inf erf allure times for each of the station. 
Tnese parameters have been computed by Regression Analysis 
for each of the stations* The exact procedure for this 
estimation of parameter is explained in appendix A. 

Table II also shews the parameter of the inter- 
failure t.ine for each of the work station. The corres- 

9 

ponding values of ’Jjf'' are also listed against them. 


TA BLE II 


W o rk S t at ion X, 

ITo , 1 

* 

*2 


Y 2 

1 

621 

615 

i 

- 630 

1.57 

2 

624 

601 

1 

" S2? 

2.96 

r- 

9 

651 

633 

1 

640 

10.96 

4 

639 

620 

1 

625 

10.71 

5 

632 

626 

1 

634 

8.67 

6 

627 

618 

1 

623 

9.23 

* 

L is 

the first r 
work st; 

: orient of 
it ion. 

the j.nterf allure 

time of the 

T a 

X 2 ' 

tlio second 
the worl 
ni.ima.te s 

moment of 
c station, 
for both 

' the Interfailure 
The values are 
X^ and Xg. 

time of 
in centi- 
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The critical value of \f 2 at 97.5 % confidence level and 

s. Ci6gr66 of freedom of 10 is 20#48« This VctluG j. $ 

2 

greater than the individual value of Y for each work 
station. This further confirms that the interfailure 
times follow a negative exponential distribution. 


5.2.2 Repair Time Distribution t 

9 

The repair time is the time taken to repair a 
failed work station. Broadly speaking, it is the time 
period after which the failed work station starts functio- 
ning again. In the present study, the repair time is the 
time for which the operator is out of the assembly line 
for his personal needs. This time is called repair time 
because after this much period of time the work station 
resumes functioning as the operator starts working. The 
time taken to satis fy the personal needs of the operators 
are nearly similar for each of them because the nature of 
their personal needs is nearly similar. Therefore it is 
assumed that the parameters of the repair time distribu- 
tion for all the work stations are equal. 

With this assumption, following are the statistics 
of the repair time distribution : 

1, Mean of the repair time = 11,10 minutes 


2 . 

3. 

4. 


Standard Deviation of the Repair time =2,85 mts. 

Y 2 = 5.75 

o / 2 

Critical value of W i,e., j cr j_^ = 


for confidence level of 97,5$ 


and degrees of freedom, of 6, 
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Cj 

Be can so \jf is less than ^ , it is concluded that the 

rer air tine follows normal distribution with a mean value 

f 

of 11.10 minutes and a standard value of 2.85 minutes. 

5.3 ALLOWANCES GI7EK TO THE OPERATORS 

Tho following allowances have been given for each 
task of tho assembly line. These values of the allowances 
are as per the current practice in the factory being 
cons. ids red* 

1* Allowance for close attention = 2% of the task time 

2, Allowance for mental strain = 1$ n 

3. Contingency allowance = 4% l! 

Total allowance =7% *' 


After alloting 
computed using 

Standard, time : 


the allowances the standard task times are 
the following expression : 

Mean value of tho observed time )i Rating 
Factor] + Allowance 


whe re , 
Rating 


Observed rate of world observed output in. the v 

factor= — S A!JJ L- 

Formal rate of Work (normal output in the . 

shift ) 


5*4 TEST LIT G OF TME METHODOLOGY OF VIF-SIF ASSEMBLY LIKE 

From the above discussions It is clear that YIF-SIF 
Assembly Lino satisfies the following assumptions : 
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1. It is a single model assembly line, 

7* Only -one operator is assigned to eacli worlc 
station, 

The task times, and lienee the work station 
times, follow normal distribution, 

4, The interfailure times follow negative 
exponential distribution, and 

5. The repair tine follow the normal distribution. 

Thus the assembly line satisfies the assumptions 

made during the development of the methodology. Therefore, 
the methodology can be applied to this assembly line. 

In Phase I, while carrying out the balancing of 
the assembly line it should be borne in mind that the 
soldering and trimming operations assoc iated with the 
components assigned to a work station have to be necessarily 
completed in the work station before the plate is passed 
onto the next station. This means that whenever a task 
of fixing a component is assigned to an operator, the task 
times for soldo ring and trimming, alongwith the task -time 
for fixing it, arc subtracted from, the time available with 
with the operator. To facilitate the understanding of the 
balancing procedure, a new term called ’component time* 
has been Introduced, Component time for component i is 
equal to the summation of component's fixing time, 
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soldering tire and the trimming time. Mathematically, 
j.t can he expressed by the following relationship : 


whe re ? 



1 



S* 


T r 

~ I 


C t 4- g t 4. rp t 

f i 1 i 

Task time for f joeing the component 

Soldering time of the i 0 * 1 component. This 
time does not include the time for bringing 
the soldering rod to the plate and the time 
for putting it back to its place. 

4*u 

Trimming tine of the i component. This 


tine does not include the time for bringing 


the pliers to the plate and the time for 
putting it back to its plate. 

The component times for the various components 
are shown in the precedence diagram on the north-east 
corner of the circles representing the various tasks. 


The methodology has been tested on the VIF-SIF assembly 
line. The results obtained are presented in the next 


chap te r • 


RESULTS AND DISCUSSIOITS 


The results obtained by using the methodology 
on trie ''ata of T VIF-SIF f Assembly Line of J .11. Electronics, 
Kanpur are presented in the following sections. In 
conformity with the generalized formulation of the problem 
given j.n chapter 3, the results are discussed separately 
for both the phases. 

6.1 PHASE I 

Table III depicts the results obtained for the 
allocation of the tasks to the various work stations by 
using the methodology for this phase. The work loads men- 
tioned against the work stations do not include the task 
time for the tasks which are identical for all the work 
stations. An examination of the work load values given in 
the table indicates that the work load is distributed 
quite uniformly among the various stations. Further, it 
yields a Balance Delay Ratio of 1.20$, a Smoothness 
Index of 13.64 and a cycle time of 725 centiminutes. 

These results are compared with the results of 
the current practice for loading the work stations. 

Table III also shows the current allocations of the tasks 
to the various work stations. It is observed that the 
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line is not balanced for these allocations because 
station 2 is much more heavily loaded as compared to the 
other work stations and therefore is a bottleneck for the 
line. From the results one notes that the assembly line 
ha s a Balance Delay Ratio of 11,30$, a Smoothness Index 
of 594.1 and the cycle time of 809 centiminutes. 

Thus it is observed that the use of the methodo- 
logy reduces the Balance Delay Ratio from 11,30$ to 

I. 20$ , indj .eating approximately 10,1$ improvement in the 
efficiencies of the work stations. The cycle time reduces 
from. 80S centiminutes to 725 centiminutes, Indicating a 

II. 5$ 'reduction. The small value of Smoothness Index 
indicates that the total idle time of the line Is distri- 
buted much more uniformly compared to the one in the 
current practice. 

The other alternative solutions representing 
equally good balances are shown in Table IV. The Smooth- 
ness Index values for these solutions range from 13,78 to 
15.56, for a Balance Delay Ratio of 1.20$. 

6.2 PHASE IX 

The output of Phase I is the Input to Phase II. 
In Phase I, the tasks are allocated to the various work 
stations while in phase II, the optimum buffer capacities 
for the buffer spaces are determined. 



u i 


" 3L :V 

lt •* i ..-/-tunr 



T ■ '• L I Vt C-.--.Tu. ) 
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Table V shows the mean and the standard devia- 
tion of the work station times. The mean of the work 
station times shown in the table also include the task 
times for the tasks which are exactly identical for all the 
work stations (e*g*, the task of bringing the soldering 
rod near the plate and putting it back to its position, 
etc.) • 

Table VI shows the values of tho number V^ 
fh 

associated with i buffer space. It may be recalled that 
V^ is the jsummation of variances of the i^ and (i+l; x 
work stations* Three groups of buffer spaces have been 
obtained using the procedure explained in section 4*2.2. 
Group 1 consists of buffer space 1, group 2 consists of 
buffer spaces 4 and 2, and group 3 consists of buffer 
spaces 3 and 5*- The three groups of buffer capacities 
yield three design variables* 

Table VII shows the results obtained after 

l 

applying the methodology for phase II. All the dejsign points 
reached during the search for optimal solution have been 
tabulated there* The starting point for the optimization 
procedure used is (0,0,0) for this table. The components 
of the design point (x-^ x 2 , x 3 ) are tho values of the 
buffer capacity sizes associated with the group 1,2 and 3 
of buffer spaces, respectively. For the results presented 
in table VII to table XI the confidence level chosen 
is 97*5$. ' Table VIII represents the results 
obtained when the methodology was applied using 
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TABUS V 


Work SJr^atJ.onJTo. 

Mean Station Time 
__(JlL_centiminute s) 

Corre sponding 
Standard Deviation 
(in centiminute) 

1 

816 

105 

2 

821 

106 

r~ 

o 

812 

131 

4 

822 

114 

5 

813 

ns 

6 

827 

138 


TABLE VI 


- v - r - 

( arranged in 
incre as ing order 


226.1 


T 


T 


271.6 ? 275.2 5 301.6 


Cor re sp o nd ing 
Buffer space 
number 


1 


0 


1 

0 

0 

0 

0 

JL 


332.0 
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TABLE VII 


SIGN OUTPUT TOTAL PERCENT UTILIZATION 
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antithetic set of random numbers. Further, the results 
obtained by taking starting design point as (24,24,24) 
are tabulated in table IX. 

Table X shows the results of the simulation runs 
when the equal buffer capacities were considered for all 
the buffer spaces. The second column of the table repre- 
sents the value of the buffer capacities. It is to be 
noted in the table that there are two rows for every value 
of the buffer capacities. The second row represents the 
results of simulation runs using antithetic set of random 
numbers. The average results of the two rows for each 
value of the buffer capacity in Table X are tabulated in 
Table XI. 

G.2.1 Discussions of the Tables : 

From Table XI, it is observed that the introduc- 
tion of buffer inventories considerably affects the rate 
as well as the cost of production from the- assembly line* 
An examination of the table indicates that the higher the 
allowed buffer capacity the greater is the rate of produc- 
tion from the assembly line. An increase of buffer 
cauacitv from aero to a value 16 for all the work 
stations (i.e.yto a total buffer capacity value of 

. , nun .'■w ' ' • 

16 x 5 = 80 for whole of the line) increases the average 
rate of production from 28.48 units per shift to 39.07 
units per shift, indicating an Increment of 26.5$. 
Further, it is observed that there is less improvement 

in the production rate after the buffer capacities for 
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all the work s J at ions have reached a value of 17, 

An examination, of Table IX indicates that 
there is an optimum value of the buffer Capacities which 
gives the minimum cost of production per unit. The cost 
of production per unit reduces from Rs. 88,20 to Rs.81.87 
if the buffer capacities for all the buffer spaces are 
increased from zero to a value sixteen. The above discu- 
ssion implies that by increasing the buffer capacities for 
all the buffer spaces from zero to a value IS, an increment 
of 26, 5i in the production rate and a reduction of 7,1$ in 
the cost of production is achieved. 

From the same table it is further observed that 
the first work station is least blocked, for there is an 
unlimited amount of units available for processing. The 
last work station is the maximum blocked and shows minimum 
utilization. This is expected because an occurrence of 
blocking at any work station increases the probability 
that it will also occur in the successiding work stations, 
Bocau.se of this cumulative effect, the last work station 
experiences the maximum blocking. 

Further it should be noted that the higher the 
buffer capacities, the greater is the utilization of the 
work stations* This is because of the fact that higher 
buffer car>acit 3 r reduces the chances of blocking and this 
in turn yields greater utilization of the work stations. 
Table VII shows the optimum solution obtained 
by using the methodology. The optimum values of 



S2 


the buffer capacities are 

1^=10, B 2 =14, B 3 =12, B 4 =14 and B 5 =12 


Corresponding to the optimum values for the 
buffer capacities, the rate of production is 3S.27 units 
per shift and the cost of production is Rs.80.72 per 
unit. These values of the production rate and the cost of 
production are the average values of the results obtained 
by u sir.” the antithetic sets of random numbers tabulated. 


in Table vn and VIII respectively. 

The simulation run of the unbalanced assembly 
line (i.e», the line before balancing in phase I) resulted 
in the following values. For the same buffer capacity 
sizes of 35 for all the buffer spaces, the rate of 
production is 32.42 units per shift and the cost of 
productj.cn per unit is Rs.B2.36, 


This shows that the methodology developed in this 


research work helps in reducing the buffer capacities 
considerably* It also results in an improvement of the 
rate of production by 17.2$ and a reduction in the cost of 
pro duct ion by 12.5$. 

It was observed from Table XI that if same values 

of the buffer capacities are used for all the buffer spaces, 

the optimum result gj.ves a production rate of 39.07 units 

Rs 

per shift and cost of production/81.87 per unit. This 
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shows that .if different values are used for the "buffer 
capacities (as is the case in the methodology) instead 
of using the sane values for all the "buffer spaces, a 
reduction of 1.5% in the cost of production and an impro- 
vement of 0.5$ in the production rate is achieved. The 
reason for this reduction of cost is attributed to the 
high buffer holding cost per unit in the last work 
station. 

The efficiency of the methodology can be improved 
by a proper choice of the initial conditions for the buffer 
capacity levels. The author has experienced that the 
initial values of the buffer sizes for a simulation run 
greatly ef '-oct the total computer time required, to reach 
the steady state of the simulation run. 


CHAPTER VII 


CONCLUSIONS AND SCOPS FOR FURTHER RESEARCH 


7.1 CCUCLTJSIOUS 

The following Inferences are drawn from this 

study : 

1. The design retliodology has been divided into two 
phases. The first phase allocates the tasks to 
various work stations so as to minimize the Balance 
Delay Ratio and Smoothness Index for the line* The 
second phase determines the optimum size of the 
buffer capacj.ties which minimizes the total cost of 
production per unit* 

2. In phase I, reallocation of tasks to the various 
work stations improves the overall efficiency of the 
line significantly. This phase gives more than one 
solutions which are equally good. Thus several 
alternative solutions are available for the line 
supervisor and he can select the one best suited to 
him from practical cons3.de rations. 

3. In phase II, introduction of buffers significantly 
affects the rate of production from the line. 
Following observations are mane in this respect : 
i) The production rate rapidly increases with the 

Increase of the buffer sizes till dhe buffer 
size Is 17. After this, the Improvement in 
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production rate due to increase in buffer sizes 
Is very ireagre* 

ii) The blocking decreases as the buffer size 
decreases, 

4. In chase xi the results of the study indicate that 
the introduction of buffers affects the per unit 
cost of assembly. With the optimum value of buffer 
obtained from the solution, the rate of production, 
shoe's an increment of 17,2% as compared to the present 
rate of production from the li ne . The improvement 
in the performance oi hie existing system is reflec~ 
ted by a reduction of the cost of production to the 
or tent of 12,5%, 

7.2 SCOPE FOR FURTHER RESEARCH 

There is ample scope for further research in 
the field of design of an assembly line system. Attempts 
may be made to relax some of the assumptions made in 
the development, of the methodology in the present work, 

A more realistic design of the assembly line would be 
achieved if the following aspects of the assembly line 
are also considered : 

1, In the methodology 5 an attempt was made to 
distribute the total work load of the 
assembly line equally among the operators, 
assuming the performance ratings of all the 
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operators to be same. However, in real 
lixe practice, the performance ratings of 
the various operators are hardly found to 
bo .the same. This feature should be 
incorporated in the methodology if the 
management is ready to pay the operators 
according to their individual performance 
ratings. This will allow the allocation of 
more work load to the operators having higher 
performance ratings. Some research work has 
been reported in this area by Mansoor |"26 j 
but it needs further investigations. 

2. The modern trends in customers 1 desire for 
variety of models and the highly competitive 
market males the management think in terms of 
mixed model assembly lines, A. generalized 
formulation of assembly line balancing 
problem should incorporate this feature. 

Some research work 37 has been reported in 
literature which deals with mixed model 
assembly line but it is very little both ■ 
in volume and substance. 

3 . In the proposed methodology it was assumed 
that all the operators are well trained and 
their performance ratings do. not change with 
time. It has been observed that In practice 
1) new operators are occasionally put on the 
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line and, 2) old operators are transferred 
to other departments. Therefore, the 
concept of 'learning 1 plays an important 
role in deciding the overall efficiency of 
the assembly line. Farther research work is 
needed to study the effect of 'learning 
process' on the behaviour of the assembly 
line . 

4. It might be interesting to carry out a 
sensitivity analysis to determine the 
dependency of the production cost and the 
rate of production upon the buffer holding 
costs and the variances of the work stations. 

Facia of the above mentioned avenues requires 
further research.. Incorporation of these features will 
help make the design of assembly line system more 
realistic . 
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APPENDIX A . 


ESTIMATION ON THE HL RAMSTER !OF AN 
EXPONENT TALLY DISTRIBUTEE FUNCTION 



The parameter of exponentially distributed set of 
emperical points is estimated by using a mathematical 
method called 'Curve Fitting** A general exponential 
relationship between the two variables x and y may be , 

♦ > ■ 

established in the Jform 

R x 
y = f 0 o 

where P Q and are parameters of the relationship bet- 
ween x and y. This relationship can also be expressed 
as 

log y = log x 

After fitting the linear curve between x and log y, 
the linear coefficient and log are computed using 
the following expressions, which are. established by 

N 


applying least-square method. 
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X” x ± y log y ± - I x ± log y 3 
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log e 
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= 1,2,. 



APP3KDIX B 


COMPUTER PROGRAM LISTING 
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% 3$C 

* PROGRAM L ISTING OF PHASE I * 

* $ 
************* ***#***#£ ***£*£ *,;,..£.£££ 


FOLLOWING IS INPUT TO THE PROGRAM 

M IS EQUAL TO TOTAL NO. OF TASKS IN THF ASSEMBLY LINE 

NOPRTR IS TOTAL NO. OF OPERATORS AVAI LABLE 

F OLOW ( I » J ) IS IMMEDIATE FOLLOWER MATRIX 

NT { I ) IS TOTAL NO. OF TASKS FOLLOWING TASK (I) 

N IS MAXIMUM OF NT(I ), $ 

NCI INITIAL CYCLE TIME TAKEN 

AVE(i) IS AVERAGE NORMAL TASK TIME FOR TASK (I) 

FOLLOWING ARE THE OUTPUT DATA 

KKK (I ) SEQUENCE OF TASKS GENERATED 

ISLAK(I) IS THE SLAK TIME LEFT AT STATION (I) 

JBASIN(I) NUMBER OF TASKS ASSIGNED TO STATION (I) 

BD(I) IS THE BALANCE DELAY FOR SEQUENCE ( I ) 

MOPR(I) IS NO. OF STATIONS NEEDED FOR SEQUENCE (I) 

OTHER IMPORTANT VARIABLES 

L FT I M { I ) TIME LEFT TO BE ASSIGNED TO THE STATION (I) 

LISTA,LISTB» FITL ST — ARE LISTS AS DISCUSSED IN THE PAPER 

WTl»WT2»WT3 f WT4,WT,REALWT ARE WEIGHTS AS STATED I NT HE PAPER 


COMMON/LABl/NSETfNT 

C0MM0N/LA82/LFTIM, LISTA,LISTB, IS LAK, KKK , NOFOL W, JR , KONTER 
C0MM0N/LAB3/WT1»WT2»WT3,WT4,WT,R£ALWT, TOTAL, AB 
C0MM0N/LAB4/F0L0W, AVE»FITLST»M»N,NQPRTR,OPR»C»NTOT 
DIMENSION JT ( 10 ',193 ) , MS£T { 100 ) » JT1 ( 100 ) , NT ( 99 ) , NSPT { IOC ) 
DIMENSION LFTIM (33 ) ,L I STA ( 5C ) , L I STB < 50) » I SLAK ( 30 ) , KKK ( 100 ) 
DIMENSION NOFQLWQOO) , JR (100 ) ,KONTER( 100 ) , WT1 ( 50 ) , WT2 ( 50 ) , WT3( 50 ) 
DIMENSION WT4( 50 ) , WT ( 50 ) , REALWT ( 50 ) »TOTAL ( 50 ) » AB ( 5f< ) 

DIMENSION SINDX(100)»MCPR(200),J BASIN (30 ) 

DIMENSION BD(2>3 f J)»NEQSX(200) , NTSK( 5 ) , NAVE (100 ) 

INTEGER FOLQW(10O,2O),AVE(1OO),FITLST(1OQ),NP{ 100 ) 

INTEGER C (60), OPR 
READ 1 , M, N» NOPT 
1 FORMAT (313) 

DO 2 I*1»M 

2 READ 3, (FQLOW( I,J) ,J=l,N) 

3 FORMAT (401 3) 

READ 9 , (NT ( I ) » I-I, M) 

PRINT 6,M 
DO T 1*1, M 

7 PRINT 8,1 , (FOLOW(!,J) f J=i,N) 
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F0RMAT(6X, I3,.IQX,?0I3) 

9 S format (?xUtotal mo, of tasks IN THE ASSEMBLY =*,I3,///,6X,*TASK*, 
11 : X » *1 MHED1 ATbLY FOLLOWING TASKS*) 

\ EAL) 4 » ( A V£ ( I ;. ) » 1 1 =1 , M ) 
read 59, MINST 
4 FORMAT ( 201 4) 

9 FORMAT (13) DC 

MINST IS THE MINIMUN STEP LENGTH THAT CAN BE TAKtM 

PRINT 44^4 

FORMAT (^X f^TASK NO • 6 X * * AVERAGE TASK TIME*) 

DO 4 A 4 5 I — = 1 , M 
PRINT 4446, 1 i , AVE( ID 
FORMAT ( 6X , 1 3 , 15X, 1 3 ) 

OF FOLLOWER MATRIX 


44 


45 

46 


M 


--DETERMINATION 
DO 5L ; I =1 , M 
DO 50 J=1 , M 
J T ( I , J ) =i) 

DO 15 0 1=1 
NK=NT ( I ) 

I F(NK.NE.O) 

GO TO 15 S)0 
CONTINUE 
KMl = i 

— J T ( I , J ) IS 


GO TO 4447 


47 


THE 


IMMEDIATELY FQLLGWIER MATRIX. IT 


J 

I 


503 


501 

502 
500 


J T ( I , J ) =1 IMPLIES THAT 
THAT J DOES NOT FOLLOW 
DO 1502 K = 1 , NK 
DO 1501 J=KM1,M 
I F { J.EQ.F OLOW ( I »K) ) GO TO 1503 

GO TO 1501 
JT ( I , J ) =i 
SO TO 1502 
CONTINUE 
CONTINUE 
CONTINUE 
PRINT 1505 


FOLLOWS I IMMEDIATELY. 
IMMEDIATELY. 


IS A BINARY MATRIX. 
JT( I , J ) =0 IMPLIES 


F ORMAT ( // , IX ,* 


FOLLOWING IS IMMEDIATELY FOLLOWER MATRIX*,//) 


506 

507 


DO 1506 I=1,M , ,, 

PRINT 1507,(JT(I,J)*' J=1 * K) 

MATRIX 

CALL TFR3MN ( M, JT ) 

,» FORMAT^ // i i X t * FOLLOWING IS TOTAL FOtLOWER 

----HERE JT(ItJ) BECOMES TOTAL FOLLOWER MATRIX 


matrix*,//) 

IT IS AGAIN 


A BINARY MATRI5 


IMPLIES THAT J FOLLOWS I AND 


511 


JT(I, J)=l 
FOLLOW I 

DO 1511 1=1, M . 

PRINT 1512, (JT( I, J ) , J=1»M) 

FORMAT (6012) TA ci<<; 

-NOFOLW(I) IS TOTAL NO. OF TASKS 


jT ( I , J ) =0 IMPLIES THAT J DOES Nl 


following task (I) 


B 


DC *.5: 3 I =,. , M 

j ; =o 

DO 1514 J = I , M 
IF! JT< I ,J ) .NE. ) Ji=J!+l 
i 52 4- CONTINUE 
1513 NQFOLW ( I ) =J1~ 

DO 111 1=1, M 

111 PRINT l«j.5,I,N0FCJLW<I ) 

1515 FORMAT (IK , *TOTAL NO. OF TASKS FOLLOWING T ASK ( *, 1 2»* ) =*» 12 ) 

READ l0.t f MSTEP,NC1 
ll?l FORMAT! 12,14) 

IP = .. 

; IP is AN INDEX TO SHOW THE CYCLE NO. BEING CONSIDERED 

C ( IP)=NCi 
li)2 CONTINUE 

PRINT 115, C{IP) 

115 FORMAT*//, IX, *CYCLE TIPS BEING CONSIDERED *,I4,//) 

NOPRTR=NQPT 
JPi = , 

* JP1 * COUNTER TO RECORD NO. OF TASKS EXCEEDING CTCLE TIME 

DO 104 I P 1 =1 , M 

IFtAVEUPI ) .GT.C(IP) ) GO TO 103 
SO TQ 1 04 

103 JP1=JP1 + L 

J I Gl=A VE ( IPX) 

J IG2=C (IP) 

J I G3=J I G1 / J I G2 
NSPT ( JPi) =JIG3+I 
M0PRTR=N0PRTR-1 

NP{ JP1)=NSPT( JP1)*C!IP )-AVE!IPl> 

NTSK(JP1)=IPI 
NAVE(IPl) =AVE(IP1) 

AVE(IP1) = C(IP)-NP( JPI) 

PRINT 116, JIG3,IP1 

116 F0RMAT(1X,I3, ♦PARA LLEL STATIONS ARE REQD. FOR TASK NO.*, 13) 

104 CONTINUE 

I F { JPI .EQ.O. AND. C ( IP ) • NE.NC1 ) GO TO 1112 

DETERMINATION OF SUM OF ALL TASK TIMES FOLLOWING EACH TASK 

DO 1015 1=1, M 

I I SQM=0 

DO 1016 J=1,M 

IF( J.EQ.I ) GO TO 1016 

I F ( JT { J, J ) .GT.D ) I ISOM=I ISOM+AVE ( J ) 

1016 CONTINUE 
1015 JR! I)=IISOM 

DO 1018 1=1, M 

1018 PRINT 101 9, I , JR ( I ) 

1019 FORMAT ( IX » *SUM OF ALL TASK TIMES FOLLOWING TASK!*, 12,*) = *, 14) 
1112 CONTINUE 

NTOT=.; 

DO 1111 1=1, M 
1111 NTOT=NTOT + AV£ ( I ) 

NCMIN=NTOT /NOPRTR+1 


IF( JPl.GT.2) GO TO 112 
I F { C ( IP)*NOPRTR.EQ.NTOT) GO TO 114 
IF(C (IP).LT.NCMIN) GO TO 112 
114 CONTINUE 

; BALANCE THE- LINE FOR THIS CYCLE TIME* CUP ) , USING COMSOAL 

CALL CCMSOL ( I P ) 

I F(OPR.GT.MOPRTR) GO TO -31 
I F(OPR.EQ.NOPRTR) GO TC .'.06 
IFIOPR.LT. MOPRTR) GO TC -08 

ID 6 IP=IP+1 1 

:(IP)=C(IP-1)-NSTEP 
IF( JP1.SQ. ' ) GO TO 102 
30 117 I <1 = 1 » JP1 
I TS=NTSK{ IK1 ) 

117 AVfc{IT$)=NAVc(ITS) 

GO TO 1, 2 

108 IP=1P+1 

NNTEP=2*N STEP 
C ( I P ) =C ( I P-1 ) -NNTEP 
IF(CUP).GT.NCMIN) GO TO 119 

1119 CONTINUE 

NNTEP=NNT£P/2 

I F (NNTEP.LT. MINST) GO TO 112 
C (IP)=C(I P)+NNT£P 
I F ( C ( I P ) . GT.NCMIN) GO TO 119 
GO TO 1119 

I F ( JP1.EQ.0) GO TO 102 
DO 118 I Kl=l » JP1 
I TS=NTSK( IK1) 

118 AVE(ITS)=NAVE(ITS) 

GO TO 102 

119 CONTINUE 

IF(JPl.EQ.J) GO TO 102 
DO 120 IK1 = 1 » JP1 
I TS=NTSK{ IK1 ) 

123 AVE(ITS)=NAVE(ITS) 

GO TO 1H2 

112 PRINT 113 

113 FORMAT ( IX* *LAST SOLUTION WAS OPTIMAL*) 

GO ,T0 132 

131 NNTEP=NSTEP/2 
IF(NNSTEP.LT. MINST) GO TO 112 
IP=IP+1 

C ( IP ) =C ( I P-l ) +NNTE P 
GO TO 102 

132 CONTINUE 
STOP 
END 



********************************************************* 

* CQMSOL IS THE MAIN SUBROUTINE WHICH BALANCES THE LINE * 

* FOR A GIVEN CYCLE TIME C(MCi) 

********* $#;}:###$##*>!!**************************** 


SUBROUTINt COMSOL(MCl) 

C OMMON/LAB^/LFTIMt LI ST A, L I STB » I SLAK»KKK , NOFOLW, JR,KONTER 
COMMON/LAB3/ WT1 1 WT2, WT3, WT4, WT , R EALWT , TOTAL, AB 

COMMON/ LAB4/ FOLOW » AVE »F I TLST » M , NjNOPRTR, OPR, C» NTOT 

DIMENSION NT(99),NSPT( 1* ) T{ . TO , cA , TCI kkkIIOP) 

LFTIMOD ),LI STA(5D) ,LISTB(50) , ISLAK(3U) ,KKK(1Q« ' 

NOFOL W( 100 ) t JRdOO ) »KONTERt 100 ) ♦ WTK 50) t WT2 ( 50 ) i WT3 ( 50 ) 

WT4 ( 50 ) , WT(50) , REALWTI 50) , TOTAL! 50) , ABt 5 -j ) 

SINDX (100 ) t MGPR(20U) ,JBASIN(30) 

BD ( 203 ) , NEQSXt 203 ) ,NTSK( 5 ) * NAVE ( 100 ) 

FOLOW ( iOOt 20 ) , AV£ ( 100) ♦ FITLST (100 ),NP(100) 

C (60) »OPR 


DIMENSION 

DIMENSION 

dimension 

dimension 

DIMENSION 
INTEGER 
I NTEGER 


MLA=5 > 

DO 3003 NUMBER=1.MLA 
SINDX (NUMBER) =10000.0 

— formation of list-a 

DO 100 K=1 » M 

I J=U 

DO 1020 1=1, M 
NK=NT ( I ) 

IF (NK.EQ.O) GO TO 1020 
DO 11? J=1 » NK 

10 I F ( FOLOW ( I , 0 ) .EQ.K ) IJ=I J+l 


020 CONTINUE 

LISTA(K) = IJ 

3PR IS THE OPERATOR 


MLEFT IS THE MO. OF 

ASSIGNMENT OF MLEFT 

MN=i 


NO. WHICH IS CURRENTLY BEING 
TASKS ASSIGNED TO OPR 
TASKS TO THE OPERATOR (OPR) 


ASSIGED TO TASKS 


QPR=1 
MLEFT =1 
I SUMM=G 

LFTIMt CPR) =C (MCI ) 
I.JS3 CONTINUE 


<=0 

FORMATION OF LIST-B 

DO 15 1=1, M 

IF(LlSTAd).EQ.O) GO TC 

30 TO 15 


16 k=k+: 

LISTB(K)=I 
15 CONTINUE 


< K=K 

FORMATION 

KK IS THE 


OF FITLIST-C 
TOTAL NO. OF 


TASKS IN LIST-B 


2103 CONTINUE 


B- 6 


<3 = 1 

DO II I K=i » KK 
K2=LISTB( IK) 

I F l AVE t K2 ) • LE .LFTI M{ OP R) ) GO TU 22 
30 TO 21 
22 FITLST(K3)=K2 
K3=K3+1 
21 CONTINUE 
<4=K3-1 

I SLAK(0PR)=C(MC1)-I$UMM 
JBASlN(OPR)=MtEFT-i 
IF{K 4.EQ.O) GO TO 2100 
30 TO 2000 
L?JO 0PR=0PR+1 

IF(OPR.LE.NOPRTR) GO T C 2102 
M OPR ( NUMB ER ) =OPR 
30 TO 3033 
Lt'2 MLEFT = 1 
I SUMM=Q 

LFTIM(0PR)=C(MC1) 


303 


GO TO 2103 

IF (K4.EQ. D KKK { MN ) =F I TL5T(K4) 
IFtK4.EQ.DG0 TO 55 

<4 is THE TOTAL NO. OF TASKS IN 

ASSIGNING WEIGHTAGES TO ELEMENTS 


FITLIST-C 
OF FITLST- 


25 


!--;^is‘TOi r NO:'o?3irKs'uWSS.G E D SO FAR 


I SUM1=0 

DO 25 IIK=1»K4 
I3=FITLST(IIK) 

I SUM1 = I SUM*.+A V£ ( 1 3 ) 

DO 26 I IK=1 »K4 

inU«>^MTI*VStKI»/FLOATIISUHl) 


I U=M-MN 
SUM2<'«0 
DO 27 1 1 K=1 1 K4 
I4=FITLST(IIK) 

A B t 14 ) =1. 0 / FLOAT ( 1 U-NOFOLW 1 1 4 ) ) 

27 SUM2=SUM2+A8tI4) 

DO 28 I IK=1 » K4 
I5=FITLST(IIK) 

28 WT2(IIK)=AB(I5)/SUM2 

I SUM3=0 

DC 29 1 1 K=1 1 K4 
I 6 =FITLST(IIK) 

29 I SUM3=ISUM3+N0F0LW ( 16 ) +a 
DO 30 IIK=IjK4 

30 WT3 (II tO* FLOAT (NOFOLW t I7)+1)/FL0AT( ISUM3 ) 

I SUM4=€ 

DO 31 IIK=!»K4 
IS=FITLST(IIK) 

31 I SUM4=ISUM4+AV£ t 18 )+JR tie J 


I9=FITL$T(IIK) 

02 WT4(I I K ) = FLOAT ( AVE ( 19) + JR CIS) )/ FLOAT ( ISUM4) 

----decision of final weight ages 

DO 40 IIK=l,K4 

43 WT( IIK)=WT1(IIK)*WT2(I IK)*WT3(IIK)*WT4( IIK) 
SUM6=Q • 0 
DD 41 IIK=1,K4 

41 SUM6=SUH6+WT(IIK) 

DO 42 I IK=l » K4 

42 REALWTIIIK)=WTUIK)/SUN6 
SUM7=O.S 

DO 43 IIK=l t K4 
SUM7=SUM7+REALWT(IIK) 

43 TOTAL(IIK)=SUM7 
X =RNOYl ( DUMMY ) 

D09C J5=lt K4 

IFIX.LE.TOTALU ) ) GO TO 51 

I F(X.LE.TOTAL(J5+l ) .AN D.X.GT. TOTAL (J5 ) ) GO TO 
i,3 CONTINUE 

51 KKK(MN)=FITLST(i) 

30 TO 53 

52 K KK ( MM ) =F ITLSTI J5+ 1 ) 

53 CONTINUE 

I F( ( J5+1) •EQ.K4 )G0 TO 54 
GO TO 55 

54 <4=J5 

GO TO 55 

55 CONTINUE 

1 1 2=KKK (MN) 

DO 70 J=l»N 

IF!FOLOW( 1 1 2 * J) • EQ *0 ) GO TO 71 
73 KONTERt J) =FOLOW( II2» J) 

71 JK1=J-1 

IF( JKl.EQ.' ) GO TO 84 
DO 72 JK2=1»M 

I F ( JK2. EQ .KKK { MN ) ) GO TO 76 
30 TO 77 
76 L I STA ( JK2 ) 

GO TO 72 
77 CONTINUE 

DO 73 J*i t JKi i ttta r it/oi- 1 

73 I F t JK2. £Q. KONTERt J) ) L ISTA ( JK2) =LISTA ( JK2 ) 1 

72 CONTINUE 
GO TO 81 

80 DO 82 JK2=ltM r 

I F ( JK2«EQ»KKK IMN) ) LISTA( JK2)-Aw 

82 CONTINUE 

81 CONTINUE 
MLEFT=MLEFT+1 
I I1=KKK(MN) 

LFTlM(CPR)=LFTIMtOPR)-AV£tIIl) 

I SUMM=ISUMM+ AVE (III) 


52 



B- 8 


C 


84 


M[, = M| ! + 1 

IF(MH.L : :.M) GO TO 3--' 

JCASIK CDPR)=PLEFT-. 

1 SLAKtCPR )=C(MC. ) — 1 SUM F 
-DOT;: PHI iAT 1 ON OF SMOOTH'JcSS INDEX 
I PP=’. 

DU ISL = :,,OPR 

IPP=IPP+I5LAK(ISL) 

I P A V = I P P/ OPR 


I SLK : 


94 


DU 94 1 SL=’„ »UPR 

I SLK=ISLK+(ISLAK(ISL)-IPAV)**2 
SMOTH= FLOAT ( ISLK) 

S I NDX( NUMBER ) =SQRT ( SMOTH ) 

PRINT 2^03, ( KKK ( I ) ,1=1 ,M) 

FORMAT (4013 ) 


PRINT 


1 I 


Vi 


1 3 


FORMAT 


( 


, SI NDXt NUMBER) 
X,*SMCOTH INDEX 


(JBASIN(I),ISLAK(I),I=1,0PR) 

*,F9«2» *TASKS ASSIGNED. . SLAK TIMES** 12t -I 


: , l 


X) ) 


93 


107 

91 

3003 


wALCULATI ON of BALANCE DfcLAY 
MOPR ( NUMBER ) = OPR 
NDUM=OPR*C(MCI; ) 
nnnn=ndum-ntot 

B D t NUMBER ) = FLOAT (NNNN ) /FLOAT ( NDUM ) 
I F ( NUMBER. EQ.i- ) GO TO 93 
D I FER=BDt NUMBER )-BD (NUMBER-1) 

I FtDIFER.LE .0.0001 ) GO TO 91 
CONTINUE 

PRINT 1: 7,60 (NUMBER) 

FORMAT (IX, *6ALANCE DELAY-*, Fes. y ) 

CONTINUE 
CONTINUE 
SMALL*!-' •' 

IF (SIND X ( NRT ) .LT.SMAL) GO TO 102 


U'2 


1 01 


[fvDX (KRT) 

) GC TO i 


LC4 


I €3 


105 


SMAL=SINDX (NRT) 

KRT =fiRT 
CONTINUE 
NCTR= 

DO X'. 3 NSUB = 1 »MLA 
DIFR=SIMDX ( NSUC )-5 

IF(ABS(DIFR) .LE.4. 

GO TO 1.3 
NCTR=NCTR+1 

NEQSX ( NCTR ) =NSUB _ , , 

I F (MCTR.EU.u ) GO TO 3 ; <-**■ 

PRX NT >, l“ 5 , SINCX < KRT > , ( l NOTR ) , NO TR-l J^Sl L OHING 

FORMAT ( / /? X , *MI N SMOOT HlfcSS INDEX-* ,FIO./, ,t rUtLuwino 


SEQUENCE NOS. 


FOR 


i ARE OPTIMAL 
3 ■'■•A CONTI MC- 
CALL ”IN(MLA,yGPR, II II ) 


THIS CYCLE TIME *,//,5«I3) 



3 - 9 


3PR=H0PR( IIII) 

print 3 : 35 , opr, iiii 

3005 FORMAT (IX, *MIN. NO. OF OPR. 
RETURN 
END 


*, 1 3, *FOR SEQUENCE *,I3) 


»**»****»...«*.....«*.«*»**«****«************************** 

* TFROMN SUBROUTINE TRANSFORMS IMMEDIATELY FOLLOWER BINARY * 

* matrty iTfT.Ji. INTO TOTAL FOLLOWER MATRIX JTlItJi 
^*^«^^*****^*****^^*^*«*****^* :it:5c3,[;i£ * ;5c;< ' : **** !lt *** :! '* 5!, * :, ‘ :,£:<C:iC ** :<C * 


10 


SUBROUTINE TFROMN? M, JT ) 

DIMENSION MSET { 100 ) , JT 1( 100 ) , JT ( 100 , 10 n ) 
DO IS J=1 ,M 
MSETt J)=0 
DO 2 J=1,M 


DO 3 1=1, M 

I F ( JT ( I , J ).NE.U) GO TO 2c 
SO TO 3 

22 MSETt J)=MSET(J)+1 

3 CONTINUE 
2 CONTINUE 

DO 4J=1»M 

IFtMSETtJ ).EQ.Q) GO TO 2 i 

GO TO 4 

23 DO 5 1=1, M 
JT ( I , J ) =0 

5 IF(I.EQ.J) JT C I , J) =1 

4 CONTINUE 
DO 6 J=l,M 

I F (MSfcT t J ) .NE.O) GO TO c. , 


GO TO 6 
24 CONTINUE 
DO T 1=1, M 
7 JTltl )=JTtI , J) 
KJ1=J~1 


I $UM=0 

DO 3 1=1, K JX 

DO 9 K=1,KJ1 

9 I SUM=ISUM+JT1 (K )*JT ( I , K) 

J T ( I , J)=I SUM 
I SUM=0 
8 CONTINUE 
JT ( J , J ) =1 
6 CONTINUE 
DO 50 I=1,M 
DO 5S J=1,M 

5t IF(JT(I,J)« GT • ) JT 1 1 , J) = 1 

RETURN 



i 4 4 4 3 »t -4 it 4? 444^4 4 4 444 4 4444444444 44 4 4 444 44 # 4 $ 4 $ # 4 ♦ ^? s$s $ 4 $ $ gfe $ 4 $ ^tjjcjjs^ 

* S UR ROUTINE' MINIMA GIVfcS THE FIRST MINIMUM ELEMENT K OF * 
- 7 HI A HR iY NOY(l). * 

4? 4 4 4 4 5 js >'c sjc 4« 4; 4 « 4 4c sfc 4 4* 4? 4« 4 $ 4 4 4 ^ .4c 4s 4 sjc 4; ^ 4c 4 4 4 s $ sfc 4 4 « # # £ # 4 * # 4 c 3 $c 4 3^ 4 4 c sfc 4 4 sgs 4c s$e 4 $ # $ 

SUBRCUT INI Mir ( N j NOY r K ) 

DIMENSION NOY ( OG) 

L , S = . C ( ■ 

00 . JK= , 0 

S F ( NCY ( JK ! .LT.LES) GO TO 2 . 

GO TO 10 ■ ' 

U S=NUY ( JK) 

K = JK 

CONTINUE 

RETURN 

&N0 



uoouu o oou oouou 


* * 

* PROGRAM LISTING FOR PHASE II * 

* * 
***.«!****************************** 

************************************************ 

* TH’' MAIN PROGRAM OPTIMIZES THE BUFFER LEVELS * 

* USING UNIVARIATE TECHNIQUE * 

************************************************ 


FOLLOWING ARE THE INPUT DATA TO THE SYSTEM 

H IS TOTAL NO. OF WORK STATIONS IN THE ASSEMBLY 

SIG(l) IS THE STANDARD DEVIATION OF (I)TH WORK -S ct! ! wnTPn 

FEXPU) I S THE PARAMETER OF NEGATIVE EXPONENTIALLY DISTRIBUTED 

INTER-FAILURE TIME FOR (DTH WORK STATION rcpair 

MUREPtl ), SIGRr.P(I) ARE THE MEAN AND STANDARD DEVIATION OF REPAIR 

TIME DISTRIBUTION FOR (DTH WORK STATION . 

Ml id) IS THF NO. OF BLFFER SPACES IN THE (I)TH GROUP OF VI 

-NJMItJ) I S T H£ { I ) T H MEMBER BUFFER SPACE OF {J)TH GROUP OF VI 

III hmm IS THC NUMBER OF ITEMS PRODUCED, AFTER WHICH OUTPUT RATE IS 

CALCULATED 

FOLLOWING ARE THE DATA RELATED TO THE COSTS OF THE SYSTEM 

a * ad | s THF fixed labour cost per shift per labour 

II CST 8 IS S THE ’labour COST OF PRODUCING AN ADDITIONAL ITEM (THIS IS 

FOR INCENTIVE PLAN) ,«• patn 

r < f.117 tc 7 nr OUTPUT HAVE OVER WHICH i HE IRLEN.iVfc IS PAID D <- T , T - 

f C3(i) IS THE COST OF HCLDING BUFFER PER UNIT, AT (DTH BUFFE 

C —Cl IS THE MATERIAL COST PER ITEM 

r QVVR IS THF FIXED OVERHEAD COST PER SHIFT 

q qVrhDI I S THE VARIABLE OVERHEAD COST PER SHIFT 

C 

c OTHER IMPORTANT VARIABLES 

i tk-aouiti t <; th k TIME OF NEXT ARRIVAL IN (DTH WORK STATION 

c IIIIItdDPR ( I ) is’ THE TIME OF THE NEXT DEPARTURE FROM (1 1™ | TATI ° N 
. t\ ic tmc t t mc at WHICH { I ) T H WORK STATION FAILS 

? y ^ AFT6R ,I,TH 

t=BF^ASS^ S&SFS STATUS 

c °SsERVin R CONm« N ?HE NC. OF UNITS PRODUCED BY (I .TH STATION TILL 

c Suin.Ji holds th-; lkits of tike tk buffer u-u has had CJ-i. 

c ITEMS IN IT SINC 2 f^I^CTmNlKF OF (UTH STATION TILL ATI, 

ci::::isT(n L ?s thI Iciuif taken^ mu.TH station to complete 


o o o o o o o o o o 


B-I2 

mm: i mi 

>1 p R \ R R f 1 ) HOLDS TH i VALUE OF THE TIME OF THE PREVIOUS ARRIVAL IN 

THL (I)TH WORK STATION 

■NPROPR ( I ) HOLDS THE VALUE OF THE TIME OF THE PREVIOUS DEPARTURE 
FROM TFT ( I ) TH WORK STATION 


FOLLOWING ARE THE DESIGN VARIABLES 

BFR(I) CONTAINS THE VALUE OF THE BUFFER CAPACITY FOR U)TH GROUP 

3 SB ( n ^ I S THE BUFFER CAPACITY BETWEEN ( I ) TH AND ( 1+4. )TH WORK STATI Of 


c* 


COM MON /OUTl /C it SOUT * SH IFT» ALAB» CST » OV VRt PRSSNt FACTQR»C3 ( 30 ) 

C OM MON / SF T / MX X , MNQ , MMM » MX 
COMMON/ FREARR/NPRARRl 30 ) ,NPRDPR( 30) 

COMMON /LAB.: /QUEfCUMQUE ,TIME,TNARV,TND PR, STATUS, ATIME 
COMHUN/LABZ/CUMUTL »CUS ERV 
C OMMON/LABB/HU, SIG 

nriywriM /CA It /PI-VP 


U ummjr*l/LA6?/MUy SI G 

C'OMMGN/FA IL/FEXP » S IGRE F, MUREP ,NFALTM, NARV rncpov f 

T NTl-GER BFR , B 38, T , AT I ME, STATUS ( 30 ) ,8(30) , CUMUTL ( 30 ) ,CUScRV ( ^ 
5ur-(io.:ci«OUE(io,30».TJME(301,TNARV(30l,TNDPR(3p) 

DIMENSION BFR (3L ) , NSQ( 20 ) ,NLL (3 r i ! , NJ1 (30, 30) , BBB { Su ) *NFALTM{3w ) 

D 1 MF NS ION (300 ) » NRT ( 3C0 ) , NARV ( 30 J , ACUMQ( 30 ) ,MU( 30 >»$IG ( 3v ) , IST(5w)» 
,T ( i 5= ? ) ,Z( 150 ) ,PCUTI L( 30) ,OUTRAT ( 100 ) , MUREP 1 30 ) »SIGRfcP ( 30) , FtXP ( 3U 

-READ THc INPUT DATA 
READ J.lG»NL,NSTEP*N 

J« j /j •* 7 ' t '1 T J 

READ XU, (NLL (I ) , I = 1 , ML) 

FORMAT ( LAI 3) 

DO 1,12 I = l»NL 

J=NLL ( I ) 

RFADnJtiMUK I»K),K = I , J) 

FORMAT (‘.AIT) 

R EADU4* SHIFT, CltSOUTA LAB, OVRHD,CST,OVVR, FACTOR 

FORMAT ( 3F.1 0.2 ) 

READ 1 15, ( C3 ( I ) t 1~ » N6 ) 

FORMAT (1AF6.2 ) 

DO 3, I =1,0 

Rt'ADEO , MUREP ( ?. ) ,SI0RLP(2),FEXP(2) 

FORMAT ( I A , TrO •_ ) 

READS 17 ,PR::SN 
FORMAT ( FI .4 ) 

READ (MU( I > , l=lf N) 

FORMAT ( 2/14) 

READ 3, (SIG( I ) * 1=1 »N) 

FORMAT ( 3: F6.i ) 

R C- AO 5 , P M H l MX , MX X , M NO 
FORMAT l LI 4) 

-CHOOSE A STARTING POINT BFR 


30) 


no 

in 

U2 

113 

1 14 
3.15 

8 

30 

UT 


SO! 


-LHUOAC * r *• 

READ 3C3, (DFR(I ),!=' 
FORMAT ( m 3 ) 


i ML) 



c FI hi) T H FUNCTION VALUE FMQ AT THE STARTING POINT 

CALL BLFL!T(NL»NLL,NJJ. r BFR,8BB) 

CALL MAIN! :ifc » ROB t FNQ ) 

PRINT i. . 1 1 FNQ * { BFR ( I ) ,I=.„,NL) 

100 FOKHATC.X ,*INITIAL OB-FK VALUE =*»F10.2»/1X»*INITIAL 
: i-->) 

mq=:. 

K = 

9 continue 


DESIGN VECTOR 


"A m 
10 


?Q3 


II 


DIRECTION TILL THE FUNCTION 


»? a 


c— 


CHOUS': THL DIRECTION OF MOVEMENT NSQ 
DO :.C I=i »ML 
I Ft I .CG.NQ) GO TO 31 
NSQ( I) = 

GO TO 1 
NSQ (I ) =/. 

CONTINUE 
N T »NSTEP 

-SEARCH IN THE CHOSEN 
CONTINUE 
DO 1'L 1=1 »NL 

BFR( I )=BFR( 1 ) +NT*NSG ( I )/ 2 
J F( BFR 1 1 ) .LE.l ) GO TO 2 lV 
GO TO n 

8 F R { I ) = BF R ( I )-NT*N5Q( I )/ 

I ZF.R0 = 1 
GO TO 13 

-FIND 1 THE FUNCTION VALUE FNT AT THE NEW POINT BFR 

CALL BUFER(NL f NLL»NJl, BFRtBBB) 

CALL MAIN (N6 » BB8» FNT ) 

PR IN" E^OO.FNT, (BFRtI) ,I=i»NL) 

FUP.'UT { ; X ? *FN”= *» F 10 , 2 * *BFR = *» 20 I 3 ) 

T F ( F‘ n • L7 • FNQ ) GO TO 12 

ART SEARCHING IN ANOFTER DIRECTION NSQ 

GO TO 13 


— S 


VALUE DECREASES 


12 J=1 

F NQ=FNT 

P R I NT ' ■ >i * FNQ » ( BFR ( I > » 1= - »NL ) 

101 FOR«ATt:X,%MHXT MIN. VALUE OF OB-FN 



1NT =*» 2 1 I A) 



.4 T = 1 ^ N T 




GO TO 1-G 



13 

I FCJ. •* G • 0 ) GO 

TO 15 


K — 1 

I Ft I ZERO- 

£Q.„ 

) GO TO ZC? 


DO 16 I -l 

»NL 

-N7*NSQ( I )/ 

16 

BFRU)=BFR(I ) 

202 

CONTINUE 



15 

NU=NQ+1 




I FtNQ.GT. 

NL) 

GO TO 17 


GO 70 9 



IT 

IF (K. "C. 

■ ) GC 

TG 25 


• =*,FI0.2»/1X»*NEXT DE»+P ~PF9 



14 


K = 


NU= ’ 

GO It; ^ 

P R I U T 2 6 t F N ( ) t ( B F R. ( I ) 1 1 = I * ML ) 
FORMAT f , X**MIN. • FUNCTION VALUE 


? I 4 ) 
STOP 
H D 


IS ' *, Flit • 3 » /1X,*DESIGN VECTOR®** 2’ 


********* 

* * 

* SUBROUTINE MAIN * 

* * 

%%%%%%%%% 

*********** 4 ******************************** 

* THIS SUBROUTINE SIMULATES THE SYSTcM FOR * 

* THE GIVEN VALUE OF BUFFER INVENTORIES * 
******* #*** 4 * 4 * ****** *********************** 


COHMON/OUT'/C l'» SCUT » SHI FT »ACAB* C ST* OVVR, PRfciSN » FACTOR *C3 130) 

C 0 M* 0 N/SFT/MXX* MNO t MMM »MX 

COMMON/PREARR/NPRARR(3C),NPRDPR(30) rT/lTltc ** xTM ~ 

COMMON/ LA 82 /QUF: ,CUMQU£ »TIME » TNARV, TNDPR, ST ATUS» ATI Mt 
CUy.'.-,'; /LAB2/CUMUTL,CUS ERV 
COMMON/SST/IST 
COMMON/ LA33/MUi 
C0MM0N/LAB5/N iir ,, _ kJ 

r nvvn s/FA I' /F c XPtSIGREF»MUREPtNFrtLTM,NARV 

DI M'ViS I UN OF ALT MI 300 ) * NFT (300 ) » NRT (300), NARV ( 30) , ACUMQ 30 
n T - v :.!!]( ( ,S IG( 3 C> * ISTl 30 ) tOUTRATC i [ -0 ) * Tt *. 5 . 0 r Z 1 * 

?L"^L"(" T 5'cUHWF(3 E "30l!llJE"!oS!?i«^L,.^TNDPR ( 3 0 . 


it SIG 


XXZZ=T 89377 
XXYY= '.85285 

y=sndy: (xxzz) 

YY=SNDY 4 ( XXYY 3 
YYY=SNCY 5 (XXYY) 


N=N6+~ 

NS® 2 *N 

FOLLOWING ARE THE 

I JKL»L 
ATI Mc=0 
JKA=t 
J KD® .1 
T (I JKL) = ‘i 
DO 1= ,0 


INIT IAL 


CONDITIONS FOR THE ASSEMBLY LINE 



STATUS (M+, )=' 

T I M i: ( I ) =0 
STATUS ( I ) = 

CUSLRVU ) = ; 

T f 1ARV ( I )=> 

N PRARR ( I )=0 
CALL NCRMAL ( I » I ST) 

TNDPRI I ) = I ST ( I ) 

NPRDPR(I>=0 
C UMUT L ( I ) = I ST ( I ) 

IFU.E C.N) GO TO 1 
QUE ( I )=0( I )/2 
IF(OUEU) ,E«.0) QUE ( I ) =1 
Ki=B( I )+i 
L = I + j, 

DO 2 J = 1 » K 1 
2 CUMQUE ( Lt J ) =0 
i CONTINUE 

PRINT 2 , (TNARV(I) ,I=1»N) 

20 F'jRMAT(lXf* NEXT ARRIVAL TIMES *,2LI6) 

PRINT 2 *. » (TNOPR(I) *1 = 1 »N) 

21 FORMAT ( LX * * NEXT DEPARTURE TIMES *,2016) 

QU5(N)=U 

NS=N-2 

N0PRT=1 

C GENERATE THE INITIAL FAILURE TIMES FOR ALL THE WORK STATIONS 

DO 81 I I I =1 »N 

CALL FAIL(III,FEXP,NFT ) 

81 NFALTM (III )=TNARV( 1 1 1 ) +NFT ( 1 1 1 ) 

100 CONTINUE 

C CHECK IF ANY OF THE WORK STATIONS FAIU 

DO $2 1 1 1 « 1 1 N 

I F ( TNARV (III) .LT.NFALT M( III )«AND.TNDPR( III ) .GE.NFALTMt III ) JG0T082 
I F(TNOFR( III ).UT.NFALTM( III ) ) GO TO 82 
I F ( QU ECU I-i. ) «EQ. ) GO ! 0 £ >C 
GO TO 43 

C GENERATE THE REPAIR TIME FDR (III)TH WORK STATION 

42 CAUL RPAIRU1 I , MUREP , $ IGREP , NRT ) 

NARV(III) =MPR«RR(I 1 1 ) + NRT till) 

TNARV (III) =?1AX0 INARV( I II), TNARV (III) ) 

TNDPRf HI ) =TNARV ( I II)-*- 1ST 1 1 1 1 ) 

STATUS (II I ) = 

GO TO 44 

42 CALL RPA1R( 1 1 1, MUR EP, S IGREP, NRT ) 

N ARV ( 1 1 1 ) =MI NC ( TNARV ( I II) , TNDPR (III)) +NRT (III) 

TNARV( III ) = MAX „> (TNARV ( III ),NARV{ III ) ) 

STATUS III 1 J 52 * 

c GENERATE THE NEXT INTER FAILURE TIME FOR (III) TH WORK STATION 

44 CALL FAIL(III,FEXP,NFT) 

NFALTM (III ) “TNARV ( III) +NFT (III) 

■82 CONTINUE 

C FIND TFE EARLIEST OCCURING EVENT 

do io 1=1, n 



\ 6 


c- 

c 


K=2*I-1 

IF(STM»lSC ). -Q.l) GO TO 3 
Z(K)=TN'\RV(I ) 
z ( K + t ) = 99 999 
GO TO 4 
2 Z (K) =99999 

Z ( K+i ) = TNOPR ( I ) 

4 com i al- 
io CONTINUE 

CALL M INC N2, Z,K9) 

C4LL SUBROUTINE -ARRIVE’ OR ’DEPART’ DEPENDING UPON WHICH bVtNT, 

ARRIVAL OR DEPARTURE, CCCURS FIRST 

DO 5 Hl=L,N 
M2=c*Hl 

5 I F ( M I • EU* K.9 ) GO TO 6 
M 3 = (K 9 +j.) /?. 

CALL ARRIVE ( MO ) 

MM 3= JKA*MMM 

I F { UUS ( N) . EQ. MM3 ) GO TC 62 

N0PRT=1 

GO TO 5 ^ 

6 CONTINUE 

CALL DEPART ( Ml, 8 ) 

MM2=JKC*MMM 

I F (QUc (M) .EQ.MM2) GO TC 61 

N0PRT=1 

GO TO 54 

61 JKD=JKC+1 

I F ( NOPRT • EQ. 1 ) GO TO 5 5 

GO TO 54 

55 PRINT 51, M1,Z (K9 ) . T -> * 

55 FORMAT! IX, ^DEPARTURE AT STATION 
NQPRT=N0PRT+I 
GO TO 8 

62 JKA=JKA+i 

I F( NOPRT. EQ. 1 > GO TO 7 
GO TO 54 

7 PRTMT 5 ■» f Z ( K S 

5„ FORKATC.X, ‘ARRIVAL ST STATION *.«.* 

NOPRT=NOPRT+i 
CONTINL " 

CONTINL ? 

MI J = I JKL*MV.M 

I F ( QUE ( U • LT • Ml J i GO TC COO 
I { I JKL+?, ) =ATI ME 


AT TIME =*,17} 


AT TIME =*»I7) 


3 

54 


11 


60 


f OR MAT ( lx 111 I M£ 1 TAKEN ' FOR PRODUCING* , 13, * SUB-ASSEM8LIES=*, 15 ) 
NOUT=GL=tN) 

•COMPUTE THE OUTPUT RA I E rl 

OUTRATU JKL +1 ) = FLO AT { NCUT ) /FLOAT \ AT iAl. i 

PRINT 6 :, OUTRATU JKL +3 ) 

FORMAT (,X,*OUTPUT RATE =*,F 10 . 4 ) 

IF(IJKL.LT.MXX) go TO 63 





64 


63 


65 


C 


97 


— t st wh'thlr ch; rate cf production has reached the steady stat 

AVOUT . : . 

A VOUT d - . 

DO 6'r JJ= ; , Mt\C 
J KK = I JKL- J.J 

A VOUT \= AVOUT l+OUTRAT ( J KK > 

A VOUT • = V OUT 2. +OUTR AT ( J KK+1 ) 

A V 0 T :„ = A V 0 U T 1 / F L C A T ( M N 0 ) 

A VQT? =A VQUT2 / FLOAT ( MNO ) 

RSHI 0 = A PS ( A VGTi-AV0T2 ) /ABSIAVDTi ) 

I FI RSHIO. GT.PRESN) GO TO 63 
GO TO 65 
I JKL=I JKL+l 

IFIlJKL.Lt_.MX) GO TO 1C J 
CONTINUE 

-TEST THE STOCHASTIC CONVERGENCE OF THE SYSTEM 
AVV=AVCTZ*SHl FT 
PRINT 66, AW 
I NDcX=0 

I F ( INDEX. EQ. 1 } GO TO 6969 
STD I =•" .■ 

DO 97 J J = 1 , M N 0 
JKi=I JKL-JJ 

STD1 = STD1 + (OUT RATI JK1 ) -AV0T2 ) )**2 
AMN=FLCAT (MNQ-i ) 


$IGMl=SQRT (STDI/AMN) 

S I GMS=S IGM3 / SQRT I AMN ) 

FIG= I S I G* TP ROB /ACC UR)* *2 
N D I F F = F NO- IFIX(FIG) 
IFINDIFF.GT.C ) GO TO 637 . 
MADD=I JKL+NDIFF 
1 JKL=I JKL + ’i 
I NDEX= 1 

I F I I JKL. LT. MACD ) GO TO i> 



6870 C Ui'iT i XL : 

66 FORMAT I // /IX , ^OUTPUT 

C COMPUTE TH~ AVERAGE 


RATE PER SHIFT = *, F8.2, // ) 

BUFF rR LEVELS AT EACH BUFFER SPACE 


00 25 NCUM=?.,N 


26 

25 

27 


a'\ 


NBl SUM=M 
N3I SM=0 

MBI =B I FCUM-1 ) +1 

DO It NSUM=i ,NBI . 

Nbl SUM=' .BI SUM+CUMQUI ( ‘-.'CUr, NSU. : I 

N B I S M = N B I SM + N SUM *C Ur OU E < * CU ^ * ij J ^ J c . 1M , , 

ACUMC I N CU H ) = FLOAT l NBIS F ) /FuuA . ( NBISUM 1-1, 

CONTINU.: „ „ 

PRINT 27, (ACUMQ(HCUM) *NCUM=2,N) 

FORMAT I IX , ^AVERAGE BUFFERS = *,liF8.3) 
-COMPUTE THE OBJECTIVE FUNCTION 
FON»v . i 
DO 30 JKE=P,N 

FON=FQN+ACUMQ(J<:) *CSUK.a-_) 


0 


mW 



B-l 8 

P I N = ._ » 7 • ’ * a V V- 7 5 . G 

C THE G • k;:ral INTERPRETATION OF PIN IS THAT IT REPRESENTS THE INCREASE 

C ? RfjDUCT IO.J PER SHIFT, CViR WHICH THE INCENTIVE MONEY IS PAID. IN THE 

C 7V ASSEMBLY LINE being CONSIDERED IT REPRESENTS THE INCREASED VALUE 

C 3F «P : R! : URMaNCQ LEVEL* P ;; R SHIFT. PERFORMANCE LEVEL EQUAL TO 75. v IS 

C CONS I LET. Li) AS THE NORMAL PACE AT J . K . ELECTRON ICS . 

I F { P I u . L L . ' .. ) GO TO 9CJ 
ACMTV=P IN*CST*FLOAT ( N ) 

GO TO 9 '2 
9 L A C N T V = ' } . ' 

902 CONTI ML : 

A.MATL = AVV^Cl 
PLAb*ALAB*FLOAT(N) 

0 VRHD‘, = FACTOR* (P LA B+ACNTV) 

F UNL = ( FON+AMATL+PL AB+O VRHDl+ACNT V+GVVR) /AVV 

PRINT 3 , FUN], 

37 FORMAT ( //iX,*FUNl=*,Fl€. 2) 

C COMPUTE THE TOTAL UTILIZATION OF EACH WORK STATION 

DO 15 1=1* N 

NT I Mt. = M AXO ( TNARV ( I ) »TNDPR{ I ) ) 

15 PCUTILt I) = ( FLOAT (CUMUTU I) ) /FLOAT ( NT I ME ) ) *100.0 
PRINT 16, (PCUT1LU ) *1 = 1, N) 

16 FORMAT ( IX, *UT I LIZAT ION =*,// ,20F8.2) 

RETURN 

END 


* * 

* SUBROUTINE FAIL * 

%%%***%***** *#%*##*#** *********************** 

* THIS SUBROUTINE GENERATES THE NEXT * 

* INTER FAILURE TIME FUR (I)TH WORK STATION * 
#£*$***#£*#))( *#£#*$$*****#***###******#******» 

SUBROUTINE FAIL ( I , FtXP ,NFT) 

DIMENSION NFT ( 3vl ) , FtXP (30) 

6 CONTINL ; 

X=RNDY5 (DUMMY ) 

R J=-ALCG( 1 »0-X ) *FEXP ( I ) 

N FT ( I ) = I F I X ( R 2 ) 

RETURN 

END 


**i ********************* 

* * 

* SUBROUTINE REPAIR * 


o o 




19 


B- 


#**#:$###*#####*#**#: #####*#**###** 

* THIS SUBROUTINE GENERATES THE * 

* REPAIR TIME FOR ( I )TH STATION * 

* * # # 4 # # # # 3); # # sjs if 4 if. # * # # * $: # # * # * # # # ❖ & =»= 


SUBROUTINE RPAIRd ,MU» SIG»NRT) 

t CONTI HU;,: 

DIMENSION NRT53GO) , MU < 30 5 , S I G ( 30 ) 

sum= : .■:: 

DU 5 JK=:,6 
R =RNDY4 ( DUMY ) 

5 SUM=SUM+R 

NRT ( I ) = IFIX(1.4I4*SIG( I ) * ( SUM-3 ... ) )+MU{ I ) 
I F { NRT { I ) . LE #' ) GO TO 6 

RETURN 

cND 


#>!( >5 #####**& #**£«*$*:* 

# * 

* SUBROUTINE MIN * 

# * 
*#*#*#*##**$!!'***#*** 

###########** s}#^**#***#*#**# ********* ###*#* 

* THIS SUBROUTINE GIVES THE MINIMUM VALUE * 

* OF THE ARRAY Z * 

a****************************************** 


SUBRCUTd MI N ( N2t Z ♦ K) 

TP THERE ARE TWO OR MORE MIN I MAS IN THE ARRAY 

WILL GIVE THE FIRST MINIMA OF THE SEQUENCE Z 
INTEGER Z(40> 

LC:S=L" '■ R 
DO :•■ JK=„»U2 
IFIZ(JK).LT.LES) GO TO 2 
GO TO 10 

2 LES=Z(JK) 

K=JK 

10 CONTINUE 
RETURN 
FND 


U 


THEN THIS SUBROUTI 


## 4 $**$# * * ##*#*#*** * ** * 

* 

SUBROUTINE NORMAL * 


¥ 

♦ 



* * 
**:}******************** 

SUBROUTINE NORMAL { L * I ST ) 

C ijMMji'l / LA BJ/MUtSIG 

D I M i ■ N S I TjM MUtiiO) »S I G( 3 C ) » 1 ST ( 30 ) 

SUM = 3 • C 
DO 3 JK =1 » 6 
RsRNCYi (DU MY ) 

5 $UM=SUM+R 

xx=(. .434*5 I G ( L > )*{ SUN-3 .*» ) 

I ST ( L ) = IF I X ( XX) -t-MU ( L) 

RETURN 
p HD 

*********************** 
* * 

* SUBROUTINE ARRIVE * 

* * 
*********************** 


*****************************^^^^*^ )!t ^! :4t ^*^* :f< ******J 

* THIS SUBROUTINE DCES ALL THc CHANGES IN THc ^STEM W I H * 
^ ad c: niic jo thF ARPIVAL OF AN ITEH IN THc { I ) T H STATION 

********************* **************** ****‘*** = ® t **** 3 * C **^ !! ^^ C *^ t * : ^ : ^ 


CONSIDERATION IS THE FIRST ONE 


SUBROUTINE ARRIVE! I ) . 

COMMON/PREARR/NPRARR( 30 |NPRDPR(30) * TTMP 

COMMON/LABT /QUE tCUMQUE *T IMc f TNARV,TNDPR» STATUS » aTI <- 

COMMOH/SST /I ST r , T ,™, 

??^nW^:j5^« ,I .f5'»^U.3»».THA R V. 3 C,,THOP R(M . 

t;.|T-G'=R ATIME, STATUS (30) 

c CHECK WHETHER THE WORK STATION UNDER 

c IN THE ASSEMBLY LINE 

r CHECK" IF" THE° BUFFER SPACE BEFORE THE STATION IS fcMPTY 

I F( QUE (I-C ) .SQ. > GO TC 26 
M=QUK I-i ) 

IOQ=M+: 

c JpSme’cumoueVatihs.oue.status.nprarr 

ifiatiIs"ge!nprbphii-li) go to 16 

CUMQUS ( I » IOQ) =CUMQUE ( I ,I0Q)+ATIM6-NPRARR * 


1 6 

17 


CUMQUti^ I » 1 00 ) =CUMQUE 1 1 ,IOQ) +ATI ME-HAXOINPRDPft 1 1-1) t NPRARR ( I ) ) 
NPRAPR ( I ) =TNARV( I ) 

ESTABLISH 3 NEXT ARRIVAL TIME AT l I ) TH WORK STATION 

STATUS ( I ) =1 
TNARV ( I ) = TNOPR( I ) 


cumqu:: n» ■ j )=cumque n»i htnarvii )-nprarr< i » 


26 


td W 



£2 m m •#£* »#* 

-WUT TILL AM ITEM IS ACDuD TO THE BUFFER SPACE BEFORE { I )TH WORK ST/ 

APR A R R ( I)=TNARV(I) 

T 1 Mf= { I ) =T'TARV (I > 

A 'i I Mil =T 1ML C I ) 

T r.i A R V ( I ) =T' 1DPR ( 1-1 )+l 
TNDPRU ) = TMOP R ( I -1 ) + IST( I) 

STATUS? I) =0 

SO TO 30 

TIME ( I )=TNARV (I ) 

-UP DATS A TIME, STATUS 
A T I M F = T I M e ( I ) 

STATUS ( I ) = ' 

ESTABLISH THE NEXT ARRIVAL TIME 
T.N AR V ( I ) =TNDPR ( I ) 

CONTINUE 

RETURN 

END 


*********************** 
* * 

* SUBROUTINE DEPART * 

* * 
***************** ****** 


****************** ********** ******************************** 

* THIS SUBROUTINE COES ALL THE CHANGES IN THE SYSTEM WHICH * 

* are DUS TO DEPARTURE OF AN ITEM FROM THE WORK STATION * 
****************** 4 ********************** ********** ********* 


SUBROUTINE DEPART ( I , B ) 

CGMMC-N/°RF.iRR/NPRARR{3C) ,NPRDPR{30 ) 

COMMON/LABl/QUE jCUMQUE ,TIME,TNARV,TNDPR, STATUS, AT I ME 
C0MMCh-/LA3’/CUMUTLtCUSERV 
C OMMOM / LA B?/ MU, SI G 
COMMON / SST/I ST 


r TNDPR(30) 


• 0MM0N/LA85/N „ 

HT EGER ATI Mt, ST AT US (SC) ,B ( 3C ) ,CUMUTL (30 ) , CUSERV ( 30 ) 

NTfcGZR QUi( 30) ,CUMQUE (30,30 ) »TIME(30 ) ,TNARV ( 30 , TND 
i iMewilON MU (2,,- ) ,SIG(20) , I ST 1 23) 

F(I.EC.N) GO TO AS 
FCQUE( n «GE. B( I)) GO TO 10 
,=QUi ( I ) 

IPDATE CUMQUE , ATIME»NPRDPR 
0Q=M+1 

— | 4 . } 

’UMOUE ( L, IOQ > =C’)MOU-: (L ,ICQ)+TNDPR( I )-MAXC (NPRDPRt I ) , NPRARR1 1 

.OMTliSL:; 

IPRDPR(I) =TNDPR(I) 


+ .U ) 



c 

c 

1 •; 

I 3 

II 

4- 

41 


15 

I'm 


7: 


IS NOh READY TO TAKE ONE ITEM FOR SERVICING 


TIME OF AN ITEM FOR THIS STATID 


TIM'. ( I ) --T NDPK { I ) 
a time=timl:( i ) 

D U £ ( I ) = 1 U I ( I ) + ■ 

-T H v. Wl.* R X S l m* ION 
STATUS! I ) = 

-= STA8LI 5H NEXT ARRIVAL AND DEPARTURE 
TNARV ( I ) = MAXG (TNARV { I ) ,TNDPR(I) ) 

CALL NORMAL l 1,1 ST) 

-UPDATE CUMUTL jCUSERV 
TNDPR ( I ) = TN AR V ( I ) + 1 ST ( I) 

CUMUTL ( I ) = CUMUTL ( I ) +IST ( I ) 

C USSR VC. ) =CUSLRV ( 1 ) +1 
GO ID 4 

I F(TNARVII) .EQ.TNDPRd )) GO TO I'L 
-SINCE N-'XT BUFFER SPACE IS FULL WAIT 
TNARV ( I )=MAXG (TNARV(l) ,TNARV( 1+1 ) ) 

GO TO 12 

TNARV ( I )= TNARV ( 1+1 ) 

TNDPR(I)=TNARV{I+J )+l 
STATUS (I ) =?-■ 

CONTINUE 

IF( STATUS 1 1 + 1 )I£q!i.ANC.TNDPR(I ) .EQ.TNDPR! 1+1 ) )TNDPR{ I+l)=TNDPR( 1+ 

y\ )-: 

CONTINUE 

RETURN 

END 


TILL AN ITEM IS PICKED UP FROM 


** ******************** 
* * 

* SUBROUTINE SUPER * 

* * 
********************** 


**********************--------t-r:r:r^::nu;r: 

* T ms SUBROUTINE GROUPS THE BUFFER SPACES IN NL GROUPS * 
•j* rc ,/ 1 c ac pypi AINSC IN SECTION 4»2»4 
**********^**********^**^ 5Se * :Jt# * ;,c ^ 5!c:<t ****^***** :!e ** #!,t * 


SUBROUTINE BUFL-RtNLt NLLt NJl » BFR, 3BB ) 
DIMENSION NLL (2*0) » NJi( 20,2-0) ,BBB(30)» BFR(30 ) 


DO 10 1 = 1 1 NL 

J=NLL ( I ) 

DO 15 Jl=2*0 


K =NJ « ( I , U i ) 
BBS (K ) = DFR( I ) 

continue: 

RETURN 

END 


o o o on 


# * # * # -r # ije & 4* # ❖ 5ft ?ft ift A 3ft 3ft sft 3ft 2ft Sft 3ft „; 


’ 5ft 3ft 3ft 3ft 3ft 3ft 3ft 3ft jft jft sft 


V 

* H 'U 3UR A f-'i LISTING FOR CHI-SQUARE TEST * 

*>'• a* 

* * * #* 4= « * $ * * # * :} 4; * * # * * # * * * jjs * * * * * -- * * * * * * * * * 


c 'GUI'S Is T H_, /.KRAY OF CaTA UNDER Ti ST 

DIN.! -SION ,'Of,S( < ) »YOBS( N ) ,XtXP(fT ) » Y-r.PCS. ) ,F{? ) , AREA (is. } 
l I ;•! 1 N Fi T * 1 ;" J YCi, J t XXUBSl i '. ) 

int.g.r xn('s,Yut‘S 

R:. A ft .vf. ,;il \I1 
FORMAT (F •. . A } 

DU 9542 J K L r" = ,54 

c ;i IS TH ■ TOTAL MUM BLR CF DATA UNDER TEST 

0 K IS THf. UTGR; ‘ OF FRC EDuM FOR THE CHI-SQUARE TOST 

R ; A.u : A, , N,K,L 
* ■•■’ FORMAT (!?,:-• 14) 

R 5 A D ..€ , UXUGSm ,1 = ) ,N ) 

>0; FORMAT lift F5.„ J 
KH?=K-4 
K P j. = K— 1 

COMPUTE THF MEAN VALUE , aVMU, AND THE STANDARD DEVIATION VALUE, 

OF THE DATA 
A MU=*. ,f. 

D U I = L , U 
i G -2 MU=M.*»»:XUBS { I ) 

- V-' J=. P 1/ FLOAT ( N) 

A S I G- • 

DO 10 I=!,H 

A S I G=ASIG+{ XXODS ( I )—AV MU ) 4*2 
S IG=SwR7 ( aSIG/FLOA T( N-i ) ) 

Df-FI TH . L “N'T. AND TFF UPPER LIMIT OF THE VALUES OF THE DATA 

DIVIDE THIS R; - Of iNTO ( K-2 ) EQUAL PARTS. ( I ) TH PART OF THIS DIV 
IS CALLED ! I ) TH GROUP 
SMALL =£ VM»I- ; . *SIG 


RIG = aVVIJ+;;. *SIG 

ST ”P= { B I G— SNA LL ) / F LOA T (K'-'I ) 

PRINT 2 , a V MU,SIG, SMALL, .%:S, STEP 

l- FOROATI.-Xf^MU^jF^.a,- SI G=*, F3 .2 , / IX , *LQWER LIMIT=*, F8.3, * 
/UPPER LIMI7=*, FO. 3 f * STEP LENGTH** , F8 . 3 ) 

C COMPUTE. PODS ID, ThF NUMBER OF DATA FALLING IN ( I )TH GROUP MADE 

DC IS I =., , K 
li. r.OKS(I)*- 
Do ' T I , ' ' 
a b=>c;cjgs ( I ) 

I F < A IS « L'.’ • SMALL ) GO TO 15 


o .. ■ j Lf «... 

16 I FtA3.GT.bIG) GO TO 11 
DO 1 ■: J=2,KMi 




E 


'• ! - ■ - S L L + F L ! r\T { J - 3 ) * 5 T £ P 

t C. =i-i:: + GY p 

i r { a * > e ■ 1 - „ » N D • A b • L T © B C 4 ) GO TO . L 

O ' TO ,„T 

.'.ijbGi J)=‘;< i r.s( j) + „ 

■ > 'J j u ; > 

CONTI -it 

Zr 7 r ; 

oiNo. < k )=F.u r ' r > ( k ) + 

cifaru . 

print , { :r)ns( ik> » ik = :,kj 

F i . K >> T { a t -<* j > U 1 '< 5 j 5 a R £ >13) 

-COM PUTS Y.-ypll)t T Ht £• XP:. CT'iD NUMBER OF DATA NALL IMG IN { 1 )TH G 
-JOB SIMPSON'S RULr TO INTEGRATE THE NORMAL DISTRIBUTION FUNCTIO 
COM PUT- Thy ;-r,p; 0 T:rD F R : GUYUCY OF DATA FALLING IN (I)TH GROUP 
-A LAM IS TH. STcP L--NGTH US’.C FOR INTEGRATION USING SIMPSON'S P.U 

a L = . . : • 

< f.=km 3 / •: 

Du b? J = , k ? 

X ATs. = C MA LL + F LG AT ( J- ,'.)*$ T C P 
YAZ= UAW.VMUJ/SIG 
XA=AB 5 ( A L ) 

Y = X L 

UN=I F I X { ( X A -X L ) / A L A M ) + 1 

f-!K 

DO f = ::,nn 

CALL FUN( YY,X) 

Y { M) =YY 

X-X+AL M 
A =0 , F 

du o i = •’ ? :k,: 

A=A+Y( 1) 

il=nk-;, 

u — 

* 1 • 

DU - ■ I = .. 5 . ‘ L ? ; 

3 =L + Y { I ) 

-ARLA( J ) I S TH OUMULATIV .. ARt A UNDF.R NORMAL PROBABILITY DISTRIB 
FONT I ON TILL LAST uRDIKATB UF U)TH GROUP 
- k- : { J 5 ~- { *L ;'•>!/?. .. ) *( Y { ; )-w ©' *A+2.T*B+Y(NN) )+&.5 

j r ! -u . . lt . *“ . j <>rfa( j > ». *:-area( ji 

C ON T INI'. 

DL : 57 J = 0 » K 2 

K = J-;„ 

k: ; ;=k/ + j 

K. f-B + J 

A P.Up S ‘<3 )= ... -\R 1 4 ( K<i- ) 

-FIJ) IS THF „^P:CTYD F :«'• vU'-KCY UF DATA FALLING IN ( J J TH GROUP 

fc;} = ur;.,h ) 

DU 5 b Js';,K1„ 

F t J ) R A (J)-'-R-A( J-1 ) 

F { K ) = , . -UU-NKN ) 

0 ! : ; . i = 5 K, 

XUXPd ) =F LOAT (N ) #F ( I ) 


r>:u T ire ),I=. ,K) 

F ' I..' " { , V , *p ( I ) , 3 -JR;-; *,;.'OF 6 .' 2 ) 

u ‘' l' 1 •' -j ( •' ■ x p ( I ) , 1=;. , k) 
n i'- n ,s are *»i 5 F 7 .i) 

^x-i 1 !!-' T H ; . 3 1 ! S:- RVttT IONS MTU CLASSES WITH AT LEAST L OBS. EXPECT 

r I ■ r 

sit. ■ l 


f, . = Y. 4 

uV r, 1 = « K 

Y,!B\ ! I ) = , 

6 Y- ;r> ( ; )= . 
j*. 

3 '! -• , K 

YfJtf.i { J ) =Y.JBS ( J) + XQ0S ( I ) 

Y . >, r’ < J ) =Y XP( j)+x::xp{ I ) 

Jh( Y Xp( J ). LT.FLQi T(U ) uO TO F 
J = J 4 

CuHTJIvU: 

c y ob : j ( j ) is th: ho. of caia in th- (J)Th group after abovesaid r 

C Y , 5'P ( J } 13 TH: NUMB? R CF EXPECTED DATA CORRESPONDING TO YOBS< J ) 

c- KNEW IS TH': R A L D:GR~fc uF FREEDOM AFTER REGROUPING 

YiJrtS < J-. J =YOBS( J-v ) +YOBS ( J ) 

Y C X P ( J - _ ) =YUXP( J-l ) +YP XP ( J ) 

KfJL W= J-', 

P R I NT 1 3 » K f-! L M t ( YQ rtS ( I ) , i = X » KNE W ) 

p K 1 ;TT 1 04 , K.U£ rt,{ Y: XP ( I ) , 1 = A , K NEW ) 

103 FORMAT (IS, TUI 3 ) 

. 4 r-R (15,3 :F6.- ) 

C COMPUTE THE CHI -SQUARE AND NO. OF DEGREE OF FREEDOM 

Crtl S !, = t, } 

00 ■) I=,,KN£W 
L L-JAT( YOBS ( I 1) 

z ,.=y:-.--;p(i > 

9 CHI S J=CUISf:+( (i-CS(Z , ,-Zt)-0.5)**2)/Z2 

N i .< = X i -1 A ~ v 

P P I ' : 7 » 5 ? Ch i i U s 00 

; .'5 F OR ’-irtT I /. ; ^ C H i ~ GUO A R r.~ $ » F\ 5 • 3 » ^ DUGRL h OF FR 'i = , - s , 15) 

P R i . 1 T ’1 ■ ■ • J 

f: FORM-' Tt , H- ) ) 

GC T ) 2 
' CuriTI <IL’ ' 

95%-S C C'NTj’ \L 
STOP 


C 

c 


T AD 

SUB ROUT i'lE FU:"’{ YY, X) 

— THIS .SUBROUTINE CaLCUL A r _S THE NORMAL 
FUNCTION V/- LU ; AT I 
Z*-,i**2 

YY- :'?{/./ . .V )/ 1.2365 
R " TUB-' 

*-rD 


OR 03-ABILITY DISTRIBUTION 


ERRATA 


Page/Line 

Printed as 

iii/8 

peo viding 

v/9 

Production Line 

8/11 

Gufahr 


25/4 

he BDR 

27/4 

where, 

is 

40/11 

V. ... 

l 

A , 



i=0 



3=0 




Read as 

p ro viding 

Production Rate 

Gut jahr 

the BDR 

where, 1 is 

as explained in 

section 4*2*2 

1 

i=l 

3=0 

t X-J. 



Call box 1 DEPART 1 


2 


To confirm this, 





Thi* book is to be returned on the 
date last stamped. 
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